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1.  Introduction 


Upconversion  processes  in  condensed  matter  are  generally  divided  into  two  types,  with  a 
distinction  being  drawn  on  the  basis  of  whether  energy  is  or  is  not  absorbed  in  the 
medium.  This  division  depends  on  whether  incident  light  initiates  population  dynamics  or 
not.  In  the  case  of  excitation  by  light  at  resonant  wavelengths  one  speaks  of 
non-parametric  upconversion,  and  relatively  efficient  redistribution  of  energy  can  occur 
within  a  medium  as  the  result  of  atomic  transitions  of  its  constituents,  generally  after  a 
delay  determined  by  the  timescale  of  internal  transitions  within  the  system.  For 
non-resonant  excitation,  parametric  upconversion  can  occur.  This  is  usually  less  efficient, 
but  the  response  is  essentially  instantaneous. 

In  this  project  we  investigated  both  types  of  upconversion  using  ultrashort  pulses,  with 
two  new  applications  in  mind.  First,  we  performed  comparative  spectroscopy  with 
continuous- wave  (cw)  and  ultrashort  pulse  lasers  to  search  for  new  channels  of  excitation 
which  might  improve  the  prospects  for  non-parametric  gain  in  the  ultraviolet  spectral 
region.  The  immense  bandwidth  of  femtosecond  pulses  is  expected  to  cause  striking 
ch^ges  in  the  multi-photon  absorption  efficiency  of  some  systems,  by  improving  the 
overall  spectral  overlap  of  upconversion  processes  involving  several  non-resonant  steps. 
Improvements  would  permit  the  realization  of  contact,  solid  state  ultraviolet  laser 
sources.  Second,  we  investigated  the  application  of  upconversion  for  improved  detection 
of  signals  in  tune-gated  imagery  of  objects  obscured  by  scattering  media.  Non-parametric 
and  parametric  upconversion  were  compared  in  this  application  for  their  ability  to  isolate 
the  wavelength  of  the  image-bearing  signal  from  scattered  light  at  the  incident  wavelength. 
Background  light  is  the  chief  limitation  to  existing  methods  of  optically  imaging  through 
highly  scattering  media. 

2.  Experimental  Methods 

2.1.  Laser  Sources 

f 

(i)  Two  existing,  mode-locked  Ti: sapphire  lasers  were  used  for  most  of  the  experiments 
reported  here.  Both  were  free-running  homemade  sources  in  which  the  mode-locking 
mechanism  was  passive  Kerr-lens  modulation.  Dispersion  compensation  with  prisms  in 
one  instrument  resulted  in  stable  pulses  of  approximately  100  fs  duration.  A  short  gain 
medium,  special  mirrors  and  third  order  conqiensation  in  the  second  instrument  resulted  in 
output  pulses  of  30  fs  duration. 

The  first  source  was  easily  inter-convertible  between  cw  and  mode-locked  operation.  As  a 
result  it  was  used  for  experiments  designed  to  compare  non-parametric  upconversion  in 
crystals  excited  by  either  continuous  or  pulsed  light.  The  second  source,  producing  30  fs 
pulses,  was  dedicated  to  imaging  experiments  requiring  the  shortest  possible  optical  gating 
signal. 
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(ii)  Finally,  an  optical  parametric  generator  and  a  white  light  continuum  generator  were 
used  as  femtosecond  tunable  light  sources  for  two  experiments  in  Pr-doped  and  Tm-doped 
materials.  These  measurements  required  two  different  wavelengths  to  initiate 
upconversion  processes. 

The  objective  in  this  part  of  the  project  was  to  test  a  new  detector  concept  requiring  two 
synchronous  pulses  of  different  wavelengths  resonant  with  transitions  of  the  chosen  ion. 
For  Pr  we  used  a  parametric  generator  pumped  by  a  regeneratively  amplified  Ti:  sapphire 
laser  system  at  800  nm.  Wavelengths  near  800  nm  and  1010  nm  were  obtained  firom  the 
pump  and  signal  beams.  For  Tm  studies,  wavelengths  near  650  nm  and  800  nm  were 
needed.  In  this  case,  the  red  beam  was  obtained  by  filtering  the  output  of  a  white  light 
continuum  obtained  by  doubling  the  pump  and  focusing  into  a  sapphire  platelet.  In  both 
experiments,  input  beams  were  then  combined  coUinearly  with  the  output  of  the  filtered 
white  light  cell  to  form  the  required  2-color  excitation  source.  Pulse  repetition  rate  of  the 
source  was  variable  between  0.1-10  MHz  and  relative  delay  between  pulse  trains  at  the 
two  Wavelengths  was  controlled  in  the  usual  way  with  a  mechanical  delay  stage. 

2.2  Time-gated  Upconversion  Imaging 

For  experiments  on  upconversion  imaging,  the  output  of  the  30  fs  Ti:  sapphire  laser  was 
divided  into  two  beams.  One  beam  was  magnified  and  used  to  illuminate  the  object.  The 
second  served  as  the  timing  reference,  or  gating  beam,  and  also  contained  a  telescope  for 
adjustable  wavefi^ont  control.  As  indicated  in  Fig.  1,  these  beams  were  combined  using 
Type  I  noncoUinear  phase-matching  to  provide  complete  spatial  isolation  of  the 
cross-correlated  signal  beam.  The  signal  beam  emerged  along  the  bisector  of  the  object 
and  gating  beams  when  the  relative  delay  of  these  beams  was  less  than  100  fs.  It  then 
passed  through  an  iris  which  blocked  light  from  second  harmonic  generation  of  the 
individual  fundamental  beams.  A  dichroic  notch  filter  provided  further  isolation  from 
scattered  light  at  the  fundamental  wavelength. 


Phosphor 
Screen 

I 

CCD  camera 

Object  Diffuser 

Fig.  1.  Imaging  apparatus  to  demonstrate  imaging  through  highly  scattering  media  using 
time-gated  upconversion  of  800  nm  light  to  a  wavelength  of  400  nm  in  a  crystal  of  LiI03. 
Delay  of  the  gating  (reference)  pulse  was  controlled  with  a  mechanical  translation  stage 
and  a  secondary  image  formed  on  a  phosphorescent  screen. 


Signal  light  was  relayed  with  an  imaging  lens  to  a  phosphorescent  screen  fabricated  on  a 
glass  slide  with  a  thin  layer  of  red  fluorescent  paint.  This  screen  converted  the  image  from 
the  SHG  wavelength  of  400  nm  to  the  600  nm  region,  thereby  permitting  efficient 
matching  to  the  spectral  response  of  the  detector  in  transmission.  This  also  retained  an 
in-line,  transmissive  geometry.  Images  were  recorded  electronically  with  a  CCD  camera 
as  a  function  of  delay  between  the  incident  beams. 

2.3  Resonance  Raman  Spectroscopy  in  Dense  Upconversion  Media 

A  cw  Tirsapphire  laser  which  was  tunable  in  the  900-1000  nm  region  was  used  to  perform 
resonance  Raman  spectroscopy  of  a  new  class  of  dense  upconversion  materials.  The  goal 
of  this  part  of  the  project  was  to  see  if  coherent  control  of  upconversion  might  be  possible, 
by  developing  an  understanding  of  the  vibrational  interactions  between  rare  earth 
impurities  and  their  nearest  neighbors  during  the  interval  of  time  immediately  following 
interaction  with  an  ultrafast  pulse.  This  information  is  contained  in  the  “local  mode” 
portion  of  the  Raman  spectrum,  and  can  be  identified  by  spectroscopy  on  and  off 
resonance  with  the  species  of  interest. 

First,  high  resolution  absorption  spectra  of  Yb3+:CsCdBr3  single  crystals  were  recorded  at 
12  K,  using  a  double  beam  spectrophotometer.  The  laser  was  then  tuned  to  wavelengths 
near  the  various  Stark  conqwnents  of  infrared  absorption  lines  of  the  Yb3+  ions.  At  each 
wavelength  the  complete  Raman  spectrum  was  recorded.  Then,  spectra  were  compared 
for  evidence  of  features  not  contained  in  the  host  crystal  Raman  spectrum,  “local  modes” 
excited  in  conjunction  with  the  cooperative  upconversion  processes  of  interest  in  this 
project.  Local  modes  were  expected  to  exhibit  resonant  enhancement  near  dipole-allowed 
absorptive  transitions  of  the  Yb  in^urity  ions,  and  indeed  one  such  mode  was  identified  in 
spectra  taken  at  liquid  helium  temperature. 

3.  Results 

3.1.  Discovery  of  Er  Upconversion  in  the  Red  Spectral  Region 

While  performing  spectroscopic  con^arisons  of  upconversion  generated  by  cw  and  pulsed 
excitation  in  Er:LiYF4  crystals,  a  new  upconversion  channel  in  the  red  spectral  region  was 
discovered.  The  fluorescence  emission  spectrum  and  the  unexpected  wavelength  which 
gave  rise  to  it  are  indicated  in  Fig.  2. 

Although  the  excitation  wavelength  corresponded  precisely  to  an  excited  state-excited 
state  transition  of  Erbium  in  a  wavelength  region  quite  isolated  from  ground  state 
absorption,  we  found  no  evidence  that  the  upconversion  mechanism  was  of  the 
“avalanche”  type.  That  is,  no  threshold  was  observed  for  upconversion  emission.  Instead 
a  smooth  quadratic  dependence  on  incident  light  intensity  was  measured  experimentally. 


Figure  2(a).  Upconversion  fluorescence  spectra  of  Er:LiYF4  excited  with  laser  light  from 
a  Tirsapphire  laser  operated  in  both  cw  and  mode-locked  conditions. 


Figxire  2(b).  Fluorescence  excitation  spectrum  of  ultraviolet  upconversion  in  Er:LiYF4 , 
showing  the  excited  state  resonance  channel  in  the  red  spectral  region. 


The  importance  of  this  discovery  is  that  it  extends  the  list  of  ions  which  can  potentially 
produce  ultraviolet  laser  emission  via  upconversion  of  red  pump  light  in  the  range  600-670 
nm,  through  avalanche  and  other  mechanisms,  to  include  Nd,  Pr,  Tm,  and  Er. 

3.2.  Demonstration  of  Time-gated  Upconversion  Imaging 

In  recent  years  a  number  of  techniques  have  emerged  which  successfully  permit  optical 
imaging  through  moderately  opaque  media.  One  technique  uses  ultrafast  time-gating  to 
select  the  least- scattered  light  from  light  pulses  transmitted  weakly  through  such  media. 
However,  an  important  limitation  of  such  approaches  arises  from  the  inescapable  presence 
of  scattered  light  at  the  fundamental  wavelength  and  intensity  imbalances.  In  an  attempt 
to  circumvent  this  limitation,  we  investigated  parametric  upconversion  of  the  image  as  a 
way  to  improve  its  contrast. 

Our  current  extension  of  the  time-gated  imaging  methodology  captures  the  entire  image  in 
a  single  shot.  No  electronic  processing  of  holographic  fringe  information  is  necessary, 
because  we  do  not  form  a  hologram  in  this  approach.  A  real  image  is  carried  on  the 
output  wave  resulting  from  cross-correlation  of  the  800  nm  reference  pulses  and  the 
object  wave  in  a  LilOj  crystal.  This  is  sinply  imaged  at  400  nm  onto  a  transmissive 
phosphor  screen  monitored  by  a  CCD  camera.  Details  of  the  experiment  and  the  results 
are  given  in  the  Appendix. 

Two  significant  advantages  have  been  identified  for  this  approach  to  optical  imaging 
through  semi-opaque  media.  The  first  is  that  high  quality  images  can  te  obtained  directly, 
without  the  need  for  any  electronic  image  processing,  through  diffusers  which  are  severe 
enough  to  eliminate  the  “ballistic”  component  in  transmitted  light.  The  second  is  that  the 
image  is  formed  as  the  result  of  an  intensity-intensity  correlation  process  which  reduces 
the  sensitivity  of  the  pattern  in  the  image  plane  to  phase  fluctuations  of  the  optical  fields. 
This  is  advantageous,  because  the  changes  in  speckle  background  which  require  averaging 
in  imagery  through  dynamic  scatterers  when  other  methods  are  used  are  not  present  with 
this  approach.  In  addition,  theoretical  considerations  suggest  that  a  version  of  our 
experiment  in  which  the  incident  reference  pulses  are  at  the  SHG  wavelength  instead  of  at 
the  fundamental  may  have  far  more  object  specificity  in  dynamic  scattering  media  (or  in 
situations  where  phase  contrast  imagery  is  desired)  than  any  existing  technique. 

3.3.  Population  Bottlenecks  in  2-Color  Upconversion 

Non-parametric  upconversion  was  also  investigated  as  a  potential  method  to  isolate  the 
detected  image  signal  from  light  at  the  fundamental  wavelength.  The  original  concept  was 
to  use  a  Pr3+-doped  crystal  or  glass  sample  as  an  ultrafast  gatable  detector  screen  in  front 
of  a  CCD  camera.  By  forming  an  image  at  the  upconversion  wavelength  instead  of  the 
fundamental  wavelength,  using  resonant  transitions  of  the  dopant  ions  to  provide 
reasonable  efficiency,  improved  image  contrast  was  expected.  In  this  case,  in  order  to 
preserve  ultrafast  gating  capability,  a  2-color  approach  was  needed,  as  indicated  in  Fig.  3. 
One  beam  tuned  to  a  ground  state  transition  wavelength  was  to  serve  as  the  object  beam 


and  the  second  beam  was  intended  to  upconvert  the  “first-excited  ions”  bearing  image 
information  fi*om  the  “first-arriving  light”  through  the  diffuser  medium. 
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Figure  3.  Schemes  for  non-parametric  upconversion  gating  with  Pr  and  Tm  ions. 

The  wavelengths  required  for  Pr:  glass  and  Pr:LiYF4  were  successfully  generated  using  the 
OPG  source  described  earlier.  It  was  intended  that  the  'G4  state  of  Pr  be  used  as  the 
intermediate  state  in  a  resonant,  two-photon  absorption  process  because  its  lifetime  was 
known  to  be  short.  Hence  population  was  not  expected  to  accumulate  in  this  state 
between  pulses  and  high  repetition  rates  could  be  anticipated. 

Unfortunately,  the  oscillator  strength  of  the  ^H4-‘G4  transition  in  both  Pr:glass  and 
Pr:LiYF4  was  so  small  that  insufficient  absorption  of  the  OPG  light  occurred.  This  was 
confirmed  in  infi-ared  absorption  scans  around  the  1  micron  region  in  both  samples  (Fig.  4 
and  5).  Blue  upconversion  emission  intensity  even  at  zero  delay  between  the  incident 
pulse  trains  was  inadequate  for  imaging.  To  solve  this  problem  we  turned  to  another  ion 
known  to  have  a  transition  wavelength  coii5)atible  with  Ti:  sapphire  excitation  which  could 
combine  with  excitation  at  a  second  wavelength  to  yield  very  bright  luminescence  firom  a 
higher-lying  state. 

The  wavelengths  for  excitation  of  Tm:LiYF4  were  generated  fi’om  the  fundamental  beam 
of  the  Ti:  sapphire  system  combined  with  filtered  continuum  light  as  described  ealier.  In 
this  case  the  excitation  dynamics  were  somewhat  different.  The  650  nm  beam  had  to  be 
maintained  at  a  low  average  intensity  to  prevent  avalanche  upconversion  from  occurring, 
and  several  intervening  states  lay  between  the  final  state  of  the  transition  and  the  initial 
state  of  the  7^  transition. 


OS 


Wftv«l«ngth  (nm) 


Figure  4.  Enlargement  of  the  unpolarized  infrared  absorption  spectrum  of  Pr:  glass  (no 
baseline  correction  applied).  No  Pr  features  are  observable  near  1.0  pm. 


Figure  5.  Enlargement  of  the  unpolarized  infrared  absorption  spectrum  of  Tm:LiYF4  (no 
baseline  correction  applied).  Spectral  features  near  1.0  pm  are  3H4-'G4  Pr  transitions. 
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In  this  case,  strong  blue  upconversion  emission  which  required  the  presence  of  both  beams 
was  indeed  observed.  However  two  new  problems  were  encountered  in  the  Tm  systenx 
First,  the  decay  3F2-3H4  was  too  slow  and  the  lifetime  of  3H4  was  too  long.  First-excited 
ions  undoubtedly  became  interleaved  in  time  with  late-excited  ions  during  the  interstate 
decay.  Second,  the  millisecond-scale  lifetime  of  the  3H4  state  caused  a  bottleneck  of 
population  bearing  “late-arriving  light”  information  from  previous  pulses  to  accumulate  in 
the  level  from  which  ultrafast  upconversion  gating  was  desired.  This  placed  a  severe 
restriction  on  the  inter-pulse  duty  cycle  which  was  acceptable.  Unfortunately  our  laser 
source  could  not  operate  below  a  repetition  rate  of  approximately  0.1  MHz.  Nor  did  the 
CCD  camera  have  enough  intensity  to  image  light  gated  from  single  pulse  events. 

Successful  implementation  of  the  concept  of  2-color,  non-parametric  upconversion 
detection  would  necessitate  that  three  basic  requirements  be  met  simultaneously.  First, 
strong  oscillator  strengths  are  needed  on  the  optical  transitions.  Second,  the  two-step 
absorption  should  ideally  not  involve  relaxation  of  the  intermediate  state.  Third,  the 
lifetime  of  the  intermediate  state  should  be  short,  so  as  to  avoid  population  bottlenecking. 
Hence  in  future  work,  different  dopant  ions  and  different  transition  combinations  need  to 
be  considered  for  increased  brightness.  Also,  host  media  which  provide  short  excited  state 
lifetimes  should  be  sought. 

3.4.  Local  Modes,  All-optical  Switching  and  Coherent  Control  of  Upconversion 
in  Dense  Media 

Non-parametric  upconversion  requires  the  interaction  of  pairs  of  dopant  ions  and  phonons 
occasionally  play  a  crucial,  controlling  role  in  some  upconversion  processes.  For  example, 
non-radiative  processes  play  a  dominant  role  in  determining  the  efficiency  of  avalanche 
upconversion  in  Tm:LiYF4[l].  All-optical  switching  has  been  demonstrated  in  systems  of 
this  type.  Hence  it  is  reasonable  to  assume  that  if  the  local  modes  of  rare  earth  ions  differ 
substantially  from  vibrational  modes  of  the  host  lattice,  some  measure  of  control  might  be 
exercised  over  upconversion  processes  by  controlling  local  mode  excitation  with  short 
pulse  techniques.  Thus,  it  should  be  possible  to  exert  coherent  control  over  upconversion 
processes,  and  it  may  be  possible  to  control  the  all-optical  switching  of  cooperative 
upconversion  using  ultrafast  pulse-shaping  or  other  techniques.  In  the  final  stage  of  this 
project  we  therefore  examined  the  local  modes  of  Yb;CsCdBr3  to  see  whether  in  this 
exemplary  upconversion  material  a  dominant  mode  could  be  identified,  and  whether  it  was 
well-isolated  from  host  vibrational  frequencies. 

The  main  results  obtained  by  Ranrnn  spectroscopy  in  YbiCsCdBrj  crystals  are  shown  in 
Fig.  6.  These  were  obtained  using  facilities  at  Los  Alamos  National  Laboratory  (details  in 
the  Appendix).  The  figure  clearly  shows  the  frequency  of  a  dominant  local  mode  of  Yb 
ions  in  this  material  to  be  very  different  from  that  of  host  modes  which  are  understandably 
more  intense.  The  Yb  absorption  spectrum  exhibits  the  same  mode.  This  suggests  that 
ultrafast  pulse-shaping  could  be  applied  to  excite  this  mode  deliberately  to  enhance 
upconversion  processes  mediated  by  phonons  in  the  first  step  of  dynamic  energy  relaxation 


in  such  an  upconversion  system.  Thus  we  have  identified  a  potential  host  material  and  a 
pathway  with  which  to  investigate  coherent  control  of  upconversion  in  dense  mpHia 


relative  frequency  [cm  '] 

Figure  6.  Raman  spectrum  of  YbrCsCdBrj  conqjared  with  the  Yb  luminescence  excitation 
spectrum  at  low  temperatures. 

Further  Raman  experiments  in  this  interesting  system  are  currently  impeded  by  lack  of 
availabHlity  of  an  infrared  Raman  spectrometer.  This  will  be  essential  for  future  work. 
Also,  experiments  to  demonstrate  all-optical  switching  in  this  material  have  been 
attempted,  but  have  not  succeeded  due  to  inadequate  cooling  design  of  available  cryostats. 
By  careful  comparison  of  techniques  with  respect  o  the  original  switching  experiments,  we 
have  determined  that  samples  must  be  cooled  by  direct  immersion  in  liquid  or  vapor 
cryogens.  Cryostats  based  on  conductive  cooling  are  inadequate  for  studies  of  this  kind  in 
which  significant  optical  absorption  causes  san:q)le  heating.  Consequently  future  work  will 
require  a  dedicated  infixed  Raman  facility  and  inqjroved  cryostats. 
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ABSTRACT 

Bistable  upconversion  luminescence  observed  in  a  "dense"  medium 
over  a  wide  temperature  range  exhibits  intrinsic  hysteresis 
governed  by  a  critical  degree  of  coherence  in  the  Lorentz  local  field. 
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Rare  earth  ions  typically  interact  weakly  with  one  another, 
even  if  they  are  nearest  neighbors  in  a  crystal.  The  inter-atomic 
coupling  strength  is  generally  insufficient  to  cause  strong 
delocalization  of  optical  excitation.  This  simply  means  that  while 
many  ions  engage  in  non-radiative  population  dynamics  as  the  result 
of  coupling  to  other  dopant  ions,  there  is  very  little  spatial 
coherence  in  their  interactions.  However,  in  every  interaction 
between  atoms  and  lasers,  considerable  optical  coherence  can  be 
generated,  and  its  dephasing  rate  can  govern  the  population  dynamics 
in  a  surprisingly  critical  way  in  dense  systems. 

Observations  in  concentrated  crystals  of  Cs3Y2Br9:10%Yb3+ 
have  revealed  intrinsically  bistable  upconversion  emission  by  Yb3+ 
ions  at  low  temperatures  [1,2].  Hysteresis  loops  of  diminishing  size 
are  observed  as  temperature  increases,  until  optical  bistability 
disappears  altogether.  A  compilation  of  intensities  at  which 
switching  occurs,  for  various  temperatures,  provides  a  diagram  (due 
to  M.  Hehlen)  of  the  boundaries  for  bistability  in  a  temperature- 
intensity  phase  space.  This  plot,  reproduced  in  Fig.  1,  provides  solid 
experimental  evidence  of  a  critical  point  for  cooperative  dynamics 
in  this  phase  space. 

The  occurrence  of  the  critical  vertex,  and  its  implications,  can 
be  understood  using  density  matrix  theory.  An  essential  ingredient 
is  the  introduction  of  migration-induced  dephasing.  The  migration 
process  is  shown  to  reduce  the  bilinear  coupling  of  population  and 
polarization  which  is  essential  to  generate  intrinsic  bistability. 
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Actual  predictions  of  cooperative  luminescent  intensity  versus 
incident  intensity  for  various  values  of  the  energy  migration  rate 
constant  p  are  shown  in  Fig.  2. 

The  curve  with  p=10  in  Fig.  2  depicts  the  critical  dephasing 
condition  at  which  the  hysteresis  evident  for  lower  p  values 
disappears.  Here,  coherence  required  to  sustain  intrinsic  bistability 
drops  to  a  critical  value.  Mathematically,  the  luminescence 
intensity  develops  a  saddle  point  as  its  two  stable  solutions  merge 
and  the  emission  curve  no  longer  folds  back  on  itself.  At  this  point 
energy  migration  amongst  excited  atoms  competes  effectively  for 
the  nearest  neighbour  population  undergoing  upconversion  dynamics, 
without  altering  excited  state  density.  Increased  migration  reduces 
the  time  available  for  adjacent  excited  atoms  to  undergo 
cooperative  upconversion.  Hysteresis  disappears  when  the  migration 
(dephasing)  rate  exceeds  the  upconversion  rate,  at  the  point  of 
"critical  coherence"  in  phase  space. 

In  dense  systems  of  interacting  atoms,  we  conclude  that  a 
bilinear  dependence  of  internal  dynamics  on  both  level  populations 
and  system  polarization  permits  optical  coherence  to  control  system 
evolution.  This  interesting  connection  between  population  dynamics 
and  coherence  in  the  local  field  opens  the  door  to  manipulation  of 
the  Lorentz  field  and  study  of  critical  behavior  in  nonlinear 
interactions  of  many  kinds  in  dense  systems. 
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Figure  Captions 

Figure  1.  A  phase  diagram  of  experimental  intensities  at  which 
switching  up  (unfilled  squares)  and  switching  down  (filled  squares) 
occurs  in  visible  cooperative  luminescence  intensity  in  Cs3Y2Br9:10% 
Yb3+  at  various  temperatures. 

Figure  2.  Calculated  pair  luminescence  intensity  (proportional  to 
the  square  of  the  population  inversion)  versus  incident  intensity  IQp 
on  resonance,  for  different  values  of  the  migration  rate  constant  p. 
The  dashed  curves  indicate  non-stationary  solutions  (unstable). 


h-lGu^ 


Dynamics 

•n  Dense  Rare  Earth  Systems 

S.C.  Rand 

UnteS  of 

(^AX.  313.647-2718;  scr@eecs.umich.edu) 

A.  Introduction 

•  pair  luminescence  and  avalanche  dynamics 

•  observations  of  Intrinsic  bistabllltv  foair  i 

critical  slowing  down  (of  avalanrho  \  ^  luminescence)  and 
9  own  (Of  avalanche  dynamics)  in  dense  media 

B.  Intrinsic  Bistability: 

•  mirrorless  hysteresis  In  10%Yh  r'e  v  o 

•  Critical  Phase  diagram  for  Yb  pair  luminLl^^e  """ 

C.  Critical  Slowing  Down; 

•  four-wave  mixing  experiments  in  1.5%Tm  livp 

3t  an  avalanche  upconversion  tranei«nn  •  optical  excitatior 

effective  excited  s^t~me  "  ''>« 

•  criticality  of  effective  decay  times  in  Tm  energy  levels 

D.  Summary; 


Cryostat 


► 


Pair 

L:»'minescence 

Sampla 


Spoctrometer 


VIS  Cooperative  Emission  Intensity 


h.  He-Ue^  cl  73  ,  0?'?^) 

Earlt^  ■  OiruJJLevs.  o.>.^A 

PRA  tl'  231^  0<?7‘»)  -* 


VIS  Cooperative  Emission  Intensity 


•  • 


9600  9800  10000  10200 

Wave  numbers  [cm 


Emission  Intensity 


— 1 - 1 - — I - f - 

200  400  600  800  1000 

NIR  Excitation  Intensity  [W  cm 


;  12-3  2-  '  ' 


72(1“0P22 


Temperature  [K] 


■^a-r  'Vt  Cs^Yir^ 


f \v^cl-CT(^  f  ?u)H  / 


Absorbance 


Effective  Response  Time  (ms) 


I  keht 


Re.(- :  (0.‘?u  *^-.1  S.  'PV.,  ,.  Rfcv.  B 

^4*  jju,t  l<  . 


Summary 


Cooparativt  Upconversion: 

•  intrinsic  bistability  occurs  in  10%Yb:Cs3Y2Br9  despite  the  weakness  of 

Yb-Yb  interactions  -  possible  role  of  selectivity  in  excitation 

•  analysis  indicates  a  balance  of  nonlinearity  and  dephasing 

•  unexpected  role  of  coherence:  rapid  dephasing  diminishes  hysteresis 

•  critical  phase  diagram  for  Yb  pair  luminescence: 

-criticality  of  the  local  field  correction- 


Avalanche  Upconversion: 

.  avalanche  dynamics  in  1 .5%Tm:LiYF4  -  critical  slowing  down  of  the 
effective  Tm  excited  state  lifetime 

•  no  evidence  so  far  of  coherence  between  atoms.  Data  explainable  with 


rate  equations  lacking  extra  degrees  of  freedom 
-criticality  of  the  local  field??-  ^ 
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Origin  of  criticality: 

•  do  these  two  closely  related  processes  exhibit  criticality  for  the 
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same  or  different  reasons? 

•  how  can  the  experimentalist  distinguish  local  field  effect  vs.  bare 

nonlinear  dynamics? 

•  coherence  vs.  incoherence? 

CAase  0  ePF^cTf  ? 
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Abstract 

We  demonstrate  image  reconstruction  of  an  object  obscured  by  a  highly  scattering  medium 
using  a  simple  time-gated  second  harmonic  generation  scheme. 


Time-gated  Harmonic  Imaging  through  Scattering  Media 
A.  Kuditcher,  M.  P.  Hehlen,  S.C.  Rand,  B.  Hoover,  and  E.  Leith 
Dept,  of  Electrical  Engineering 
1311  Beal  Avenue,  University  of  Michigan 
Ann  Arbor,  MI  48 109- 1 1 1 2 
Tel:  313-763-6810 
Fax: 313-647-2718 

A  significant  limitation  to  existing  optical  techniques  for  imaging  through 
obscuring,  highly  scattering  media  is  the  contrast  limitation  due  to  non-image-bearing, 
scattered  light  at  the  frequency  of  incident  beams.  For  example,  in  the  holographic  time¬ 
gated  technique  [1],  such  light  contributes  to  the  average  background  intensity,  thereby 
decreasing  contrast.  In  the  present  work  we  upconvert  an  image  using  time-gated,  second 
harmonic  generation  of  incident  light  from  an  ultrashort  pulse  laser  to  reduce  the 
interference  from  non-image-bearing  scattered  light  at  the  input  frequency.  We  extend 
earlier  results  [2]  by  capturing  an  entire  image  without  any  scanning  and  consider 
implications  of  this  technique  for  imaging  through  dynamic  scatterering  media. 

The  output  of  a  Ti:sapphire  laser  emitting  mode-locked  pulses  of  about  30  fs 
duration  was  divided  into  two  beams.  One  beam  was  magnified  5x  and  collimated  with  a 
Newtonian  telescope  to  serve  as  the  object  beam.  In  this  way,  the  object  beam  was  assured 
of  being  smaller  in  the  focal  region  than  the  reference,  thereby  improving  spatial  resolution 
for  imaging  applications  while  simultaneously  making  the  beam  large  enough  for  the 
insertion  of  coarse-grain  objects.  A  variable  delay  line  was  placed  in  the  second  beam  for 
timing  control.  A  1  mm  diameter  wire  was  placed  in  the  first  beam  at  the  object  plane  and 
two  ground  glass  screens  inserted  at  positions  approximately  1  cm  and  10  cm  beyond  it  to 
obscure  the  object  completely.  The  second  beam  was  similarly  collimated  and  aligned  to  be 
parallel  to  the  first.  Both  beams  were  focused  by  a  single  30  cm  lens  into  an  oriented 
crystal  of  LiI03  at  an  external  crossing  angle  of  approximately  4  degrees.  The  apparatus  is 
shown  in  Fig.l. 

When  the  crystal  was  rotated  to  the  Type  I  phase  matching  angle  for  800  nm,  both 
fundamental  input  beams  simultaneously  showed  efficient  harmonic  conversion  irrespective 
of  the  reference  delay.  This  was  somewhat  unexpected  because  the  inter-arm  angle 
exceeded  the  calculated  angular  acceptance  at  800  nm  for  monochromatic  beams,  but  did 
not  interfere  with  signal  detection  since  the  spatial  separation  was  more  than  adequate  to 
isolate  any  desired  beam.  When  the  delay  line  was  adjusted  for  zero  inter- arm  delay,  an 

additional,  cross-correlated  signal  beam  emerged  along  the  bisector  of  the  other  two.  The 
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intensity  of  this  signal  depended  critically  on  delay  and  bore  upconverted,  time-gated  image 
information  as  the  result  of  the  intensity-intensity  correlation  process  in  the  nonlinear 
crystal.  An  imaging  system  was  used  to  relay  images  fomied  on  a  phosphorescent  screen 
to  a  CCD  camera  for  data  acquisition. 

The  cross-correlation  of  the  scattered  object  wave  and  the  reference  beam  with  both 
diffusers  in  place  is  given  in  Fig.  2(a).  Its  full  width  at  half  maximum  was  53  fs.  The 
cross  correlation  without  the  diffusers  shown  in  Fig.  2(b)  was  only  slightly  wider  (57  fs), 
suggesting  that  pulse  stretching  attributable  to  the  scatterers  is  small.  However,  the 
wavevectors  of  severely  delayed  light  lie  outside  the  acceptance  cone  for  cross  correlation 
in  our  method,  and  fail  to  contribute  to  this  signal.  With  the  scatterers  in  place,  the  object 
could  not  be  imaged  from  light  in  the  object  beam  in  any  plane  at  the  fundamental 
frequency.  The  time-gated  harmonic  signal  however  yielded  the  results  shown  in  Fig.  3 
without  any  electronic  image  processing.  The  results  clearly  show  a  coarse-grain  image  of 
the  wire  located  behind  two  widely  separated  ground  glass  plates. 

A  simple  phosphorescent  screen  fabricated  by  spraying  a  thin  layer  of  red 
fluorescent  paint  on  a  glass  slide  sufficed  as  a  transmissive  screen  which  converted  the 
image  at  400  nm  to  the  red  spectral  region.  This  resulted  in  a  good  match  to  the  response 
curve  of  the  CCD  camera  while  retaining  the  signal-to-noise  advantage  of  working  at  a 
wavelength  different  from  that  of  the  fundamental  beams.  Image  field  data  were  recorded 
at  upconverted  wavelengths  for  various  delays.  Highest  fidelity  was  attained  at  the  shortest 
delays.  Rapid  degradation  of  image  quality  occurred  as  the  gate  delay  was  increased  from 
0  to  50  fs. 

In  summary  we  have  demonstrated  that  time-gated  upconversion  can  be  used  in  a 
non-scanning  geometry  for  background- free  imaging  of  stationary  objects  with  millimeter 
resolution.  Our  results  were  obtained  using  a  crossed  beam  configuration,  but  collinear 
(inline)  geometries  are  feasible  in  Type  II  geometries.  No  image  processing  was  necessary 
to  obtain  the  results  presented  here.  Images  were  obvious  to  the  naked  eye  at  short  gate 
delays.  This  scheme  is  expected  to  be  well-adapted  to  imaging  problems  involving  moving 
scattering  media  because  the  nonlinear  process  utilized  to  reconstruct  the  image  is  an 
intensity  correlation  which  should  be  insensitive  to  phase  variations  of  scattered  fields. 
References^ 
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Figure  Captions: 

Figure  1.  Apparatus  for  imaging  through  highly  scattering  media  using  time-gated,  second 
harmonic  generation.  Double-headed  arrows  indicate  lenses. 

Figure  2.  Cross-correlation  of  the  object  wave  and  the  reference  beam  (a)  without  any 
diffusers  in  place  to  show  the  pulsewidth  prior  to  "stretching",  and  (b)  with  two  diffusers 
in  place  to  obscure  the  object. 

Figure  3.  Upconvened  image  of  the  object  plane  (a)  without  the  object  present,  and  (b) 
with  a  1  mm  diameter  wire  in  the  object  plane,  obscured  by  two  frosted  glass  screens  and 
recorded  in  transmission  with  a  CCD  camera  from  a  secondary  image  on  a  semi-transparent 
red  phosphorescent  screen  sensitive  to4(X)  nm. 
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Figure  1. 
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abstract 


Pronounced  electron-phonon  coupling  is  observed  for  the  ^F7/2^^F5/2  4f  transitions  of 
doped  into  CsCdBr3.  A  comparison  of  the  Raman  spectrum  and  the  luminescence 
excitation  sideband  accompanying  the  ^F7/2(0)->^F5/2(2’)  crystal-field  transition  reveals 
vibrational  properties  of  the  [YbBre]  coordination  unit  that  differ  markedly  from  those  of 
the  CsCdBrj  host.  In  particular,  the  vibronic  transition  associated  with  the  totally 
symmetric  [YbBr6]  stretching  mode  appears  as  a  very  weak  feature  at  191  cm"'  in  the 
Raman  spectrum,  whereas  the  totally  symmetric  [CdBrs]  stretching  mode  of  the  CsCdBr3 
bulk,  which  appears  as  a  strong  feature  at  162.5  cm'*  in  the  Raman  spectrum,  is  only 
weakly  discernible  in  the  sideband.  This  is  direct  evidence  for  a  large  contribution  from 
{YbBre]  local  modes  and  a  small  contribution  from  bulk  modes  to  the  vibronic  intensity. 
The  intensity  of  the  local  mode  is  enhanced  by  approximately  a  factor  of  two  in  Raman 
spectrum  when  the  laser  is  tuned  into  resonance  with  the  ^F7/2(0)->^F5/2(2’)  absorption  of 
Yb^^,  providing  direct  confirmation  of  its  assignment.  The  observation  of  the  first  and 
second  members  of  a  Franck-Condon  progression  for  both  the  local  and  the  bulk  indicates 
that  a  A  process,  rather  than  a  M  process,  induces  the  vibronic  intensity.  Huang-Rhys 
factors  of  5iocai=0.010±0.002  and  5buik=0.15±0.03  were  determined  from  the  data,  and 
reflect  quite  different  electron-phonon  coupling  strengths.  These  results  suggest  that 
multiphonon  relaxation  of  excited  electronic  states  proceeds  by  the  excitation  of  local 
modes  of  [YbBr^]  followed  by  energy  transfer  to  bulk  modes  of  the  lattice,  possibly 
through  a  nonlinear  coupling  mechanism  which  is  discussed  briefly. 


I.  INTRODUCTION 


The  optical  spectra  of  non-centrosymmetric  solids  doped  with  trivalent  rare-earth 
(RE^"^)  ions  typically  are  dominated  by  electric-dipole  induced  electronic  transitions 
between  the  states  of  the  4f  shell.  These  electronic  transitions  can  be  accompanied  by 
vibronic  sidebands  which  are  a  result  of  interactions  of  the  4f  electronic  system  with  the 
various  vibrational  modes  of  the  surroundings  of  the  ion.  With  the  exception  of 
Pm^^,  vibronic  sidebands  have  been  observed  for  all  the  ions  in  a  variety  of  solids 
[1-4].  The  integrated  intensity  of  the  vibronic  sideband  relative  to  its  zero-phonon  line  is 
usually  low,  although  vibronic  transition  probabilities  approaching  or  even  exceeding 
zero-phonon  transition  probabilities  have  been  found  for  certain  transitions  in  some 
compounds  [5-8].  Recently,  Gd^^  and  Pr^^-doped  crystals  received  particular  attention 
for,  besides  having  advantageous  energy  levels,  they  represent  the  extrema  of  the 
unexplained  minimum  of  the  electron-phonon-coupling  strength  in  the  middle  of  the 
lanthanide  series  [9]. 

Vibronic  sidebands,  although  usually  weak  and  often  ignored,  are  just  one  of  the 
many  manifestations  of  electron-phonon  coupling  in  rare-earth-ion  doped  solids. 
Multiphonon  relaxation  of  excited  states  and  phonon-assistance  in  non-resonant  energy 
transfer  are  other  ubiquitous  phenomena  of  key  importance  to  excited-state  dynamics  and, 
ultimately,  to  the  role  rare-earth  ions  play  in  their  many  applications.  Knowledge  of  the 
electron-phonon  coupling  strength  as  well  as  the  spatial  extension  of  the  electron-phonon 
interaction  is  therefore  a  prerequisite  to  a  detailed  understanding  of  the  energy  flow  in 
electronic  dynamics  involving  the  vibrational  system.  Direct  information  on  the  spatial 
extension  of  electron-phonon  interactions  in  particular  is  difficult  to  obtain,  since  in  most 
luminescent  materials  such  as  RE^*-doped  YAG,  YLiF4,  LaCb,  or  Cs3Lu2Br9,  in  which 
the  RE^^  ion  substitutes  for  a  chemically  similar  ion,  the  vibrational  properties  of  the 
nearest  surroundings  of  the  RE^^  ion  are  almost  identical  to  those  of  the  host  lattice.  One 
exception  is  the  class  of  compounds  containing  well-defined  molecular  species.  The 
ternary  chloride  (NH4)2GdCl5  for  example  contains  NH4^  molecular  units  in  the  second 


coordination  sphere  of  Gd^"^,  and  electron-phonon  coupling  beyond  the  first  [GdCl^] 
coordination  sphere  to  vibrational  modes  was  observed  in  the  ^P7/2->*S 

luminescence  sideband  [10].  In  (NH4)3YCl6;Eu^^,  electron-phonon  coupling  to  those  high 
frequency  NH4^  vibrational  modes  in  the  second  coordination  sphere  was  suggested  to 
account  for  the  unusually  low  Eu^"^  luminescence  quantum  yield  of  this  material  [10]. 
Rare-earth-ion  doped  compounds  such  as  CaF2,  CdF2,  CsMgCh,  or  CsCdBrs,  in  which 
ions  substitute  for  divalent  metal  ions,  constitute  another  class  of  materials  for 
which  electron-phonon  coupling  beyond  the  first  coordination  sphere  may  be  directly 
observed.  Various  defect  structures  are  formed  in  these  materials  as  a  result  of  charge 
compensation.  These  defects  typically  have  different  symmetry  properties  which,  in 
addition  to  the  different  mass  and  charge  of  the  ion,  may  significantly  alter  the 
vibrational  properties  compared  to  the  respective  host  lattice  site.  In  this  context,  the 
vibrational  system  may  be  regarded  as  being  composed  of  the  local  vibrational  modes  of 
the  RE^^  defect  embedded  in  the  bulk  modes  of  the  host  lattice. 

In  this  paper  we  present  a  detailed  study  of  vibronic  coupling  in  Yb^^-doped 
CsCdBrj.  The  comparison  of  the  Raman  spectrum  with  the  luminescence  excitation 
spectrum  of  the  sideband  accompanying  the  ^F7/2(0)^^F5/2(2’)  crystal-field  transition  of 
Yb^^  reveals  strong  coupling  to  the  vibrational  modes  of  the  first  Yb^^  coordination 
sphere  and  weak  contributions  from  coupling  to  modes  of  the  second  coordination  sphere 
and  beyond.  In  particular,  resonant  Raman  spectra  are  used  to  assign  a  weak  feature  in  the 
Raman  spectrum  to  the  local  mode  associated  with  the  highest-energy  vibration  of  the 
first  Yb^^  coordination  sphere.  The  appearance  of  two-phonon  replica  of  both  the  local 
and  the  bulk  mode  as  well  as  of  the  respective  combination  mode  in  the  vibronic  sideband 
is  direct  evidence  for  a  Franck-Condon-type  process  inducing  the  vibronic  intensity.  The 
results  provide  insight  into  the  spatial  extension  of  vibronic  interactions  and  to  the 
various  steps  required  for  the  nonradiative  relaxation  of  electronic  excited  states  of  rare- 
earth  ions  in  solids. 
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II.  EXPERIMENT 


Optical  quality  crystals  of  Yb^^’-doped  CsCdBr3  were  grown  using  the  Bridgman 
technique  from  99.5%  cationic  purity  (c.p.)  CsBr  (dried  in  vacuum  at  250  °C)  and  CdBr: 
(obtained  from  99.998%  c.p.  CdO  and  HBr,  and  sublimed  at  440  °C),  and  YbBra.  The 
anhydrous  YbBrs  was  obtained  in  a  four-step  process  by  formation  of  (NH4)3YbBr6  from 
evaporating  a  47%  HBr  solution  of  Yb203  (99.999%  c.p.)  and  NH4Br  (99.8%  c.p.),  then 
drying  of  (NH4)3YbBr6  in  a  nitrogen  stream  at  200  ®C,  decomposition  of  (NH4)3YbBr6  to 
NH4Br  and  YbBr3  at  420  °C  in  vacuum  [11],  and  final  sublimation  of  YbBr3  at  1040  °C 
in  a  boron-nitride  ampoule.  The  CsCdBr3:Yb^'^  sample  was  obtained  from  cutting  a 
crystal  along  planes  containing  the  crystallographic  c  axis  and  applying  an  optical  polish 
to  the  surfaces.  The  total  ytterbium  concentration  of  the  sample  used  for  this  study  was 
1.47±0.04  mol%  as  determined  by  inductively  coupled  plasma  optical  emission 
spectroscopy  (ICP-OES).  Since  Yb^^  is  easily  reduced  to  Yb^^  and  moreover,  Yb 
substitutes  for  Cd^^  when  doped  into  CsCdBr3,  typically  both  Yb^"^  and  Yb^^  will  be 
present  in  a  Yb-doped  CsCdBra  crystal  grown  by  the  method  described  above.  The  yellow 
color  of  the  crystal  is  indication  for  the  presence  of  some  Yb^^,  and  therefore  the  Yb^^ 
concentration  is  likely  to  be  smaller  than  the  above  value. 

Polarized  absorption  spectra  were  recorded  on  a  double-beam  spectrophotometer 
with  the  sample  cooled  to  19  K  by  a  closed-cycle  helium  refngerator.  An  argon-ion  laser 
pumped  Ti:  sapphire  laser  was  used  as  excitation  source  for  the  near-inirared 
luminescence  and  luminescence  excitation  spectra.  The  sample  luminescence  was 
dispersed  by  a  1-m  single  monochromator  and  detected  by  a  cooled  InGaAs  diode  using  a 
lock-in  amplifier.  The  unpolarized  Raman  spectra  were  obtained  using  the  polarization- 
scrambled  488.0  or  514.5-nm  line  of  an  argon-ion  laser.  The  Raman  spectrum  was 
independent  on  the  choice  of  either  one  of  these  wavelenghts.  The  scattered  light  was 
collected  in  a  backscattering  arrangement,  dispersed  by  a  1-m  double  monochromator, 
and  detected  by  a  photomultiplier  and  a  photon  counter.  Unpolarized  resonant  Raman 
experiments  using  an  argon-ion-laser  pumped  Ti:sapphire  laser  as  excitation  source  were 


recorded  in  a  similar  geometry  with  the  sample  cooled  to  12  K  by  a  closed-cycle  helium 
refrigerator.  The  scattered  light  was  dispersed  by  a  triple  monochromator  and  detected  by 
a  liquid-nitrogen-cooled  CCD  detector. 

Except  for  the  absorption  and  the  resonant  Raman  spectra,  all  experiments  were 
carried  out  at  either  78.5  or  15.5  K  with  the  sample  mounted  in  a  cold-finger  cryostat. 


III.  RESULTS  AND  DISCUSSION 
A.  ^t^Energy  Levels  of  Yb^  in  CsCdBrj 

The  diamagnetic  salt  CsCdBrs  is  a  member  of  the  inverse  perovskite  family  with 
the  hexagonal  CsNiCl3  structure  [12].  It  crystallizes  in  the  space  group  and  consists 
of  linear  chains  of  face-sharing  [CdBr6]'‘'  units  with  the  chains  being  ananged  along  the 
crystallographic  c-axis  and  the  Cs^  ions  occupying  high-symmetry  sites  between  the 
chains.  It  was  first  shown  by  EPR  measurements  that  trivalent  rare-earth  ions  (RE^"^) 
preferentially  are  incorporated  as  charge  compensated  RE^'^-vacancy-RE^^  ion  pairs 

t 

[13,14]  having  an  RE^'"-R£^'^  intra  ion-pair  distance  of  approximately  6.0  A  [14].  Several 
other  RE^^  minority  sites  have  since  been  identified.  The  point  syriimetry  for  a  RE^^  ion 
on  a  Cd^*  site  is  Csv  [15-19]  thus,  for  odd  4/electron  systems,  the  multiplets  are 
split  into  (2  J  + 1)/2  Kramers  doublets. 

In  its  [Xe]4/^  electron  configuration,  only  has  two  multiplets,  i.e.  a  'F7/2 
ground  state  and  a  excited  state  which,  in  Csv  symmetry,  are  split  into  4  and  3 
crystal-field  levels  labeled  (0),  (1),  (2),  (3)  and  (O’),  (!’),  (2’),  respectively.  The  three 
prominent  lines  observed  at  10122,  10138,  and  10599  cm'*  in  the  19-K  absorption 
spectrum  (Figure  1)  can  consequently  be  assigned  to  the  (0)— ►(0’,r,2’)  zero-phonon 
transitions  of  Yb^^,  respectively.  Each  of  the  transitions  is  accompamed  by  a  multi-peak 


sideband  which  results  from  electron-phonon  interactions  discussed  in  Sections  III.  B  and 
C.  The  (0)->(0’,r,2’)  energies  are  very  similar  to  those  measured  in  CssYbiBrg  (10119, 
10146,  10590  cm  ,  respectively  [20])  in  which  has  a  6-fold  bromide  coordination 
geometry  closely  resembling  the  one  of  in  CsCdBr3:Yb^^.  From  the  luminescence 
spectrum  (Figure  2),  the  (O’)^(0, 1,2,3)  transitions  are  found  at  10122,  10007,  9955,  and 
9658  cm’’,  along  with  the  (r)->(0, 1,2,3)  transitions  which  are  observed  as  a  result  of 
thermal  population  of  (T)  at  15.5  K.  These  energies  indicate  a  ground-state  splitting  of  0, 
1 15,  167,  and  464  cm  ’  which  again  is  quite  similar  to  that  found  in  Cs3Yb2Br9  (0,  1 14, 
140, 441  cm’’  [20]).  Both  the  ground  and  excited-state  splittings  reflect  the  dominance  of 
the  octahedral  over  the  trigonal  component  of  the  crystal  field,  since  the  (0)-(l,2)-(3)  and 
the  (0’,r)-(2’)  splittings  are  much  larger  than  the  (l)-(2)  and  (O’)-(r)  splittings.  The 
and  ^F5/2  crystal-field  energies  determined  here  agree  with  thoM  reported  in  an  earlier 
study  [21]  except  for  (2)  and  (2’).  A  small,  unexplained  shift  of  the  (0)->(2’)  transition  to 
higher  frequency  with  increasing  temperature  was  observed.  Therefore,  the  slight 
discrepancy  in  these  energies  between  the  two  studies  most  likely  result  from  the  different 
temperatures  (10  K  vs.  15.5  K)  at  which  the  spectra  were  recorded. 


B.  Bulk  and  Local  Modes  in  CsCdBrscYb^ 

t 

The  Raman-active  first-order  normal  modes  (corresponding  to  phonons  near  AK) 
around  the  Brillouin  zone  center)  of  the  CsCdBrs  host  can  be  derived  from  a  nuclear  site 
group  analysis  of  the  unit  cell  [22,23].  In  the  space  group  of  CsCdBrs  with  two 
formula  units  in  the  primitive  unit  cell,  the  Cs,  Cd,  and  Br  atoms  occupy  sites  of  » 
and  symmetry,  respectively  [24].  These  three  sites  give  rise  to  the  modes 
A2u+Big+Eiu+E2g,  A2u+B2u+Eiu+E2u»  and  Aig+A2g+A2u+Big+Biu+B2u'*'Eig+ 
2Eiu+2E2g+E2u,  respectively,  with  the  modes  A2g+2Big+Biu+2B2u+2E2u  being  silent  and 
A2u+Eiu  being  the  acoustic  modes  [23].  From  the  remaining  active  optical  modes 


2A2u+3Eiu+Aig+Eig+3E2g,  which  we  refer  to  as  bulk  modes,  Aig+E|g+3E2g  are  Raman 
active.  All  these  modes  are  observed  in  the  unpolarized  Raman  spectrum  shown  in  Figure 
3.  Both  the  energies  and  relative  intensities  of  the  Raman  transitions  are  in  agreement 
with  previously  reported  Raman  spectra  of  undoped  CsCdBr3  [25, 26]. 

Pronounced  vibronic  sidebands  are  observed  in  both  the  absorption  spectrum 
(Figure  1)  and  the  luminescence  spectrum  (Figure  2).  In  the  discussion  below,  we  focus 
on  the  (0)-^(2’)  zero-phonon  line  and  its  vibronic  sideband  since  it  is  well  separated  in 
frequency  from  all  the  other  transitions.  The  luminescence  excitation  spectrum  of  the 
vibronic  sideband  accompanying  this  (0)-^(2’)  zero-phonon  line  (Figure  3)  reveals  quite 
a  complex  structure,  with  many  more  lines  spanning  a  wider  frequency  range  than  in  the 
Raman  spectrum.  Such  a  pronounced  qualitative  difference  between  the  Raman  spectrum 
and  the  vibronic  sideband  is  typically  not  observed  for  rare-earth-ion  doped  solids  in 
which  the  RE^^  ions  substitute  for  chemically  alike  ions.  We  do  not  attempt  to  analyze  all 
the  features  of  the  sideband  shown  in  Figure  3.  Rather  we  focus  on  the  highest-energy 
transitions  which  are  distinct  and  well  isolated.  In  the  isostructural  compounds  CsCoCh 
and  CsMgClj,  the  Aig  bulk  mode  was  found  to  have  the  highest  frequency  of  all  the 
Raman  and  IR-active  optical  modes  [27].  Also,  no  modes  with  frequency  higher  than  the 
Aig  bulk  mode  were  found  in  an  analysis  of  vibronic  sidebands  in  CsCdBrsiCo^"^  [25].  We 
therefore  assume  that  the  Aig  Raman  transition  at  162.5  cm"'  represents  the  highest- 
energy,  first-order  bulk  mode  in  CsCdBrs.  The  highest-energy,  first-order  feature  in  the 
vibronic  sideband  on  the  other  hand  is  observed  at  191  cm'*,  a  frequency  close  to  the 
value  (190  cm*')  which  was  observed  for  the  highest-energy  mode  in  both  Cs3Yb2Br9  and 
Cs3Er2Br9  crystals  with  six-fold  bromide  coordination  [20,28].  As  shown  in  the  inset  of 
Figure  3,  the  191  cm''  mode  in  the  vibronic  sideband  is  also  discernible  as  a  very  weak 
feature  in  the  Raman  spectrum.  Since  there  are  no  second-order  Raman  transitions  which 
could  give  rise  to  a  transition  at  this  frequency,  we  assign  this  weak  Raman  transition  to 
the  Aig  local  mode  associated  with  the  [YbBr6]  coordination  unit.  This  assignment  is 
directly  confirmed  by  resonant  Raman  experiments.  Figure  4  shows  the  Raman  intensity 
of  the  191-cm''  local  mode  {I  i)  relative  to  that  of  the  162.5-cm''  bulk  mode  (/g)  as  a 


function  of  excitation  frequency  and  compares  I l! I b  the  unpolarized  (0)^(2’) 
absorption  spectrum.  The  increase  oi  I  ij  I  g  by  approximately  a  factor  of  two  on  the 
•  (0)-»(2’)  resonance  is  direct  evidence  that  the  191-cm‘‘  mode  in  the  Raman  spectra  is  due 

to  a  local  mode  associated  with  Yb^"^. 

The  relative  frequencies  of  the  Aig  local  and  the  Aig  bulk  modes  can  be  explained 
^  on  the  basis  of  a  simple  model.  These  modes  represent  totally-symmetric  stretching 

vibrations  of  the  M-Br  bond  (M=Cd,Yb),  with  frequencies  proportional  in  a  first 
approximation  to  where  k  is  the  force  constant  and  /r  is  the  reduced  mass.  Their 

relative  frequencies  should  therefore  be  given  by 


^YbBr  _  M^YbBr/A:dBr 
^CdBr  V  ^CdBr/^T)Br 


Assuming  the  six  M-Br  bonds  to  be  equivalent,  each  of  the  Br'  ions  shares  of  the 
charge  of  the  metal  ion  in  the  [MBr6]  coordination  unit,  i.e.  and  of  an  elementary 
charge  for  [YbBre]  and  [CdBr6],  respectively.  Further  assuming  the  force  constant  to  be  a 
linear  function  of  the  ionic  bond  strength,  a  ratio  of  ^vbBr/^CdBr  =  1-132  is  predicted 
from  Eq.(l).  This  value  is  within  4%  of  the  experimental  ratio  191/162.5=  1.175, 
providing  further  support  for  the  mode  assignments. 

The  fact  that  the  Aig  local-mode  frequency  is  so  similar  to  the  respective 
frequency  in  Cs3Yb2Br9  but  significantly  different  from  the  A|g  bulk-mode  frequency 
indicates  a  high  degree  of  flexibility  in  the  JcdBrf  j  linear  chains  of  CsCdBra. 
Obviously,  forces  exerted  by  the  host  lattice  do  not  constrain  the  [YbBr6]  normal-mode 
frequencies  to  be  the  same  as  CsCdBrs  bulk-mode  firequencies.  The  Yb^"^  defect  site  is 
apparently  able  to  accommodate  local  forces  by  bending  of  adjacent  metal-halogen-metal 
bond  angles.  This  flexibility  of  the  JcdBrf  j  linear  chains  is  reflected  not  only  in  .the 
large  angular  distribution  of  the  crystallographic  c-axis  of  1112®  measured  in 
CsCdBrsrPr^^  using  hole-burning  Stark  spectroscopy  [29]  but  also  in  the  fact  that  metal 
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ions  over  a  wide  range  of  charge,  size,  and  chemical  characteristics  can  be  doped  in 
significant  amounts  into  crystals  of  the  CsCdBr}  family  [24]. 

In  Figure  3,  a  weak  transition  in  the  excitation  sideband  at  164  cm'*  coincides 
with  the  strong  Aig  bulk  mode  observed  at  162.5  cm'*  in  the  Raman  spectrum.  This  weak 
transition  may  be  either  of  a  vibronic  or  an  electronic  nature.  Cooperative  electronic 
transitions  on  the  Yb^"^  ion  pairs  in  this  material  are  likely  and,  for  example,  two  excited 
Yb^""  ions  cooperatively  relaxing  from  their  excited  state  in  a  radiative 

[^F5/2(i’),^F5/2(j’)]->[^F7/2(m),^F7/2(n)]  transition  can  give  rise  to  intense  green  cooperative 
luminescence  in  CsCdBrsiYb^^  [21].  In  analogy,  cooperative  electronic  transitions  of  the 
type.  [^F7/2(i),^F7/2G)]->[^F5/2(ni’),^F7/2(n)]  can  occur  in  the  near  infrared  spectral  region 
[20].  From  the  energy  levels  calculated  in  Section  III.A,  the  transition 
[^F7/2(0),^F7/2(0)]-^[^F5/2(2’),^F7/2(2)]  is  cxpectcd  at  10766  cm*',  i.e.  at  a  frequency  167 
cm'*,  above  the  (0)->(2’)  zero-phonon  line.  We  do  not  have  enough  experimental 
information  to  assess  the  contribution  of  this  electronic  process  to  the  transition  intensity 
around  164  cm  *.  However,  it  is  likely  to  be  small,  since  the  observation  of  a  progression 
in  a  vibrational  mode  of  approximately  164  cm'*  (see  Section  III.C)  indicates  significant 
vibronic  character  of  the  164  cm'*  transition.  Moreover,  a  transition  at  162  cm  *  was  also 
observed  in  CsCdBr3:Pr^'^,  which  has  a  different  ground-state  splitting  [30].  Therefore, 
this  transition  arises  mainly  from  coupling  of  the  (0)-^(2’)  electronic  transition  either  to 
the  Aig  bulk  mode  in  the  second  coordination  sphere  or  to  an  odd-parity  local  mode  of  the 
[YbBr6]  coordination  unit  with  a  frequency  similar  to  the  Aig  bulk  mode. 

In  the  vibronic  sideband  accompanying  the  ^H4(0)-»^Po(0)  zero-phonbn  line  in 
CsCdBr3:Pr^^,  the  Ai,  local  mode  of  the  [PrBr6]  coordination  unit  appears  at  194  cm'*, 
and  another,  fairly  broad  vibronic  transition  appears  at  162  cm'*  [30].  The  slight  increase 
of  the  Aig  local  mode  firequency  from  191  cm'*  in  CsCdBrsrYb^"^  to  194  cm  *  in 
CsCdBr3:Pr^*  is  expected,  due  to  the  smaller  mass  of  Pr^"^  relative  to  Yb^^.  However  the 
fact  that  the  164  cm  *  mode  frequency  in  CsCdBr3:Yb^'^  does  not  increase  by  going  to 
CsCdBr3:Pr^'^  is  indication  for  this  transition  not  being  a  result  of  vibronic  coupling  to  a 


local  mode  but  arising  from  vibronic  coupling  of  the  (0)->(2’)  electronic  transition 
beyond  the  first  [YbBr6]  coordination  sphere  to  the  Aig  bulk  mode. 

These  results  clearly  show  that  vibronic  coupling  occurs  almost  exclusively  to  the 
local  modes  of  the  first  coordination  sphere  with  minor  contributions  from  the  second  or 
higher  coordination  spheres,  i.e.  a  maximum  vibronic  interaction  radius  of  ^3  A.  In  terms 
of  electron-phonon  interactions  enabling  non-radiative  relaxation  processes  such  as 
multiphonon  relaxation  or  phonon-assisted  energy  transfer,  this  conclusion  implies  that 
most  of  the  energy  released  by  these  processes  is  first  transferred  to  the  local  modes  of  the 
first  coordination  sphere.  The  subsequent  distribution  of  this  vibrational  energy  among 
the  various  bulk  modes  requires  resonances  of  those  bulk  modes  with  the  local  modes,  a 
condition  which  may  not  be  well  satisfied  for  defect  materials  such  as  RE^^-doped  CaF:, 
CdF2,  CsMgCl3,  or  CsCdBr3.  Time-resolved  Raman  studies,  such  as  those  carried  out  for 
Cr‘*'‘-doped  forsterite  [31],  could  provide  further  insight  into  the  details  of  the  non- 
radiative  relaxation  dynamics  in  these  materials.  They  might  also  answer  the  question  as 
to  whether  or  not  the  lack  of  any  strong  resonances  between  local  and  bulk  modes 
contributes,  together  with  the  overall  low  frequencies  of  phonons  in  CsCdBr3,  to 
suppression  of  non-radiative  relaxation  of  RE^"^  excited  states  in  this  material. 


C.  Electron-Phonon  Interaction 

Two  mechanisms,  usually  referred  to  as  M  and  A  process  [32],  can  mediate  the 
interaction  of  the  electronic  and  the  vibrational  system.  The  M  process  induces  vibronic 
intensity  through  admixture  of  even-parity  wavefunctions  {e.g.  from  the  4f'‘5d 
configuration)  to  the  odd-parity  4f  wavefunctions  by  coupling  to  an  odd-parity  (IR- 
active)  vibration  [3,  33-35].  Such  vibrationally-induced,  forced  electric-dipole  transitions 
lead,  for  example,  to  vibronic  transitions  with  a  frequency  differing  by  exactly  one 
vibrational  quantum  from  the  zero-phonon  line.  The  A  or  Franck-Condon  process  on  the 


other  hand  induces  vibronic  intensity  through  slightly  different  equilibrium  structures  of 
electronic  ground  and  excited  states  [36].  In  contrast  to  the  M  process,  the  shift  of  the 
excited-state  relative  to  the  ground-state  potential  surface  along  one  of  the  configurational 
coordinates,  measured  in  terms  of  the  Huang-Rhys  factor  S,  gives  rise  to  phonon  replica 
in  the  respective  mode.  In  the  harmonic  approximation,  the  probability  of  a  vibronic 
transition  from  an  electronic  state  a  and  vibrational  state  n  to  an  electronic  state  b  and 
vibrational  state  m  is  given  by 


Pa„-*b„  =  Pab\(xbM\Xain)f 


(2) 


where  P^i,  is  the  purely  electronic  transition  probability,  which  is  identical  for  all  the 
vibrational  states  m  and  n,  and  |(/lfi(w)|ylfa('*))|  the  vibrational  overlap  integral  or 
Franck-Condon  factor,  which  is  a  measure  of  the  electron-phonon  coupling  strength  [36]. 
At  zero  Kelvin,  the  relative  transition  probability  for  the  /th  phonon  replica,  given  by 
Pao-^b,  =  Pab\{xb(0\xai^)f  [37] 


Pab 


S‘e~^ 


(3) 


Numerous  studies  in  recent  years  have  found  that  earlier  asstimptipns  that  the  Huang- 
Rhys  factor  is  zero,  implying  that  the  M  process  dominates  the  RE^^-ion  vibronic 
intensities,  were  often  not  justified-  The  observation  of  vibronic  transitions  in  the  two- 
photon  spectra  of  SrMo04:Pr^^  [38]  and  of  two-phonon  replica  in  Na5La(M04)4:Pr^^ 
(M=Mo,W)  [39]  provides  direct  evidence  for  the  important  role  that  the  A  process  can 
play. 

As  shown  in  Figure  5,  the  bulk  and  local  modes  vb  and  vt  observed  at  164  and 
191  cm  '  in  the  vibronic  sideband  of  the  luminescence  excitation  spectrum  appear  as  very 
weak  two-phonon  replicas  at  323  and  383  cm'',  respectively.  That  is,  they  appear  at 
approximately  twice  thier  respective  ftmdamental  frequencies.  This  is  direct  evidence  for 


a  A  process  inducing  the  vibronic  intensity.  This  result  is  consistent  with  selection  rule 
considerations  since  both  the  modes  vb  and  vi  were  assigned  to  Aig  modes,  which  as 
even-parity  modes  cannot  enable  M  processes.  The  two-phonon  frequencies  are  at  «  1.96 
and  «  2.01  times  the  one-phonon  frequencies  for  the  modes  vb  and  vl,  respectively, 
indicating  the  highly  harmonic  nature  of  the  excited-state  potential  surface.  Note  that  the 
combination  mode  vg  +  is  also  observed. 

By  using  Eq.(3)  it  is  possible,  in  principle,  to  calculate  the  Huang-Rhys  factors  S 
from  the  ratios  of  transition  probabilities  Pa^^bJ^ab  and  Pa^^b^  j  Pab  derived  from 
experimental  transition  intensities  /  in  the  luminescence  excitation  spectrum  (Figure  5). 
Under  the  experimental  conditions  used  for  measuring  the  luminescence  excitation 
spectra,  however,  it  is  likely  that  the  strongly  absorbing  zero-phonon  line  is  partially 
saturated,  i.e.  its  intensity  is. not  proportional  to  Pab  •  Consequently,  the  above  two  ratios 
would  be  overestimated,  a  complication  in  numerous  past  studies  [4]. 

Instead,  we  assume  that  the  weakly  absorbing  vibronic  transitions  are  unaffected 
by  saturation  {i.e.  I  <x  P),  and  we  only  calculate  the  ratio  I by  using 

Eq.(3)  and  the  spectra  of  Figure  5.  For  the  modes  vg  and  vi  we  then  obtain 
S Bulk  =0.15±0.03  and  =0.010±0.002,  respectively.  Note  that  Sguik  >  ^Locd  • 
result  reflects  the  fact  that  the  relative  intensity  of  the  modes  vb  and  vi  in  the  one-phonon 
and  two-phonon  replica  (Figure  5)  is  reversed.  Although  both  Huang-Rhys  factors  found 
here  are  within  the  typical  range  of  values  reported  for  other  RE^'^-doped  solids  [7,39], 
they  indicate  quite  different  electron-phonon  coupling  strengths  of  the  two  modes. 

Whereas  the  configurational  coordinate  for  the  Aig  local  mode  is  most  likely  the 
(average)  Yb-Br  bond  length,  the  configurational  coordinate  describing  the  coupling  to 
the  bulk  mode  in  the  second  coordination  sphere  or  beyond  is  not  directly  obvious.  The 
relatively  large  value  of  S indicates  the  occurrence  of  local  structural  changes  upon 
excitation  of  the  respective  vibronic  transitions  and,  on  the  basis  of  both  the  one¬ 
dimensional  natiue  of  CsCdBra  and  the  structure  of  the  dominant  Yb^^  pair  center  in 
CsCdBr3,  we  can  speculate  on  the  nature  of  this  structural  rearrangement.  From  our 
results  on  the  spatial  extension  of  electron-phonon  interactions  (Section  III.B)  we 
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conclude  that  vibronic  coupling  is  not  only  confined  to  within  the  linear  chains  in 
CsCdBr3  but  more  specifically  does  not  extend  beyond  the  vibrational  modes  of  the 
nearest-neighbor  coordination  units  in  the  chain.  With  the  dominant  Yb^"^  pair  center 
being  a  Yb^^-vacancy-Yb^^  arrangement,  the  Yb^"^  ion  therefore  typically  has  a  [CdBr^] 
coordination  unit  on  one  side  and  a  [Vacancy-Br^]  coordination  unit  on  the  other  side 
along  the  linear  chain.  Since  the  vacancy  is  likely  to  be  easily  compressed,  the  vibronic 
excitation  of  the  Aig  bulk  mode,  which  involves  stretching  of  the  Cd-Br  bonds,  may  push 
the  Yb^"^  ion  towards  the  vacancy  and  thus  induce  changes  in  bond  lengths  and  angles  of 
the  bromide  ions  linking  the  Yb^"^  and  the  Cd^"^  ion.  The  overall  result  of  this  vibronic 
excitation  would  be  the  introduction  of  coupling  between  local  and  bulk  modes  and  a 
reduction  of  the  Yb^'^-Yb^'^  distance.  While  there  is  currently  no  experimental  evidence 
supporting  this  picture,  we  wish  to  point  out  that  EPR  measurements  on  a  sample  strongly 
excited  into  the  bulk-mode  vibronic  transition  at  162.5  cm  *  could  provide  direct 
information  on  the  magnitude  of  this  structural  rearrangement  [14]. 


IV.  CONCLUSIONS 


The  incorporation  of  into  CsCdBrs  leads  to  the  formation  of  [YbBre] 
coordination  units,  the  vibrational  properties  of  which  differ  markedly  from  those  of  the 
host  lattice.  From  a  comparison  of  Raman  spectra  which  probe  host  vibrations,  and  Yb^"^ 
luminescence  excitation  spectra  which  provide  a  probe  of  local  vibrations,  we  have 
shown  that  vibronic  coupling  occurs  predominantly  to  modes  of  the  first  coordination 
sphere,  and  that  vibronic  contributions  from  the  second  or  higher  coordination  sphere  are 
weak.  The  spatial  extent  of  electron-phonon  interactions  is  estimated  to  be  approximately 
3  A.  These  results  all  have  important  consequences  for  phonon-assisted  energy  transfer 
and  multiphonon  relaxation  of  excited  states.  The  present  research  has  shown  that 
vibrational  energy  release  during  relaxation  of  the  electronic  excited  state  of  Yb^^  occurs 
through  excitation  first  of  local  modes,  and  subsequently  via  energy  transfer  through 
resonances  with  bulk  modes  which  complete  the  nonradiative  relaxation  process.  For 
materials  such  as  RE^'^-doped  CsCdBr3,  where  local  and  bulk  mode  frequencies  differ 
markedly,  the  significant  detuning  of  many  such  potential  resonances  undoubtedly 
suppresses  some  relaxation  pathways. 
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FIGURE  CAPTIONS 


Figure  1:  Unpolarized  ^F7/2->'F5/2  absorption  spectrum  of  CsCdBr3:Yb^"'  at  19  K.  The 
labels  (0)  and  (0’,r,2’)  refer  to  the  crystal-field  levels  of  the  ^F7/2  ground 
state  and  ^F5/2  excited-state  multiplet,  resp)ectively. 

Figure  2:  Unpolarized  ^F5/2->‘F7/2  luminescence  spectrum  of  CsCdBr3:Yb^'^  at  15.5  K. 

The  labels  (0,1, 2, 3)  and  (0’,!’)  refer  to  the  crystal-field  levels  of  the  ^F7/2 
ground  state  and  ^F5/2  excited-state  multiplet,  respectively. 

Figure  3:  Unpolarized  Raman  spectrum  of  CsCdBr3:Yb^^  with  the  incident .  laser 
wavelength  at  514.5  nm  (top  trace)  and  unpolarized  luminescence  excitation 
spectra  obtained  by  monitoring  the  (0’)->(2)  transition  (Figure  2)  at  1004.5 
nm  (bottom  traces).  For  the  Raman  and  excitation  spectra  the  frequency  is 
normalized  to  the  incident  laser  and  the  (0)-^(2’)  zero-phonon  line  (Figure 
1),  respectively.  The  assignments  of  the  Raman  transitions  are  taken  from 
Ref.  [25,  26].  The  inset  shows  an  enlargement  of  the  Raman  spectrum  in  the 
1 80-200  cm*'  range. 

Figure  4:  Relative  unpolarized  Raman  intensity  of  the  bulk  mode  at  162.5  cm  '  (/g) 
and  the  local  mode  at  191  cm*'  (7^)  as  a  fimction  of  the  incident  laser 
frequency  in  CsCdBr3:Yb^'^  at  12  K.  The  enhancement  of  the  ratio  /7g  on 

the  (0)— K2’)  resonance  is  approximately  a  factor  of  two.  Part  of  the 
unpolarized  absorption  spectrum  (from  Figure  1)  is  shown  for  comparison. 

Figure  5:  Unpolarized  luminescence  excitation  spectra  of  CsCdBrsiYb^*.  vg  and 

denote  the  bulk-mode  and  local-mode  frequency  respectively,  and  the  left  and 
right  insets  show  the  one-phonon  and  two-phonon  replica  of  both  these 
modes,  respectively. 
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Two-beam  coupling  measurements  of  an  avalanche 
that,  despite  the  resonant  nature  of  the  excited-5tat( 
dispersive.  This  surprising  characteristic  is  shown  to 
off-resonant  process  dominating  the  resonant  response. 

In  recent  years  simple  implementations  of  two-bedm 
coupling  interactions^  have  been  introduced  for  pre¬ 
cision  characterization  of  nonlinear  coefficients^  and 
for  studying  slow  radiative  and  nonradiative  decay 
processes  within  luminescent  centers  in  solids.^  Here 
we  report  results  of  nondegenerate  two-beam  coupling 
<TBC)  mediated  by  delocalized,  cooperative  dynamics 
of  impurity  atoms  on  avalanche  upconversion  transi¬ 
tions^  and  use  them  to  investigate  an  earlier  indication 
from  research  by  Ni®  that  avalanche  nonlinearities  are 
primarily  dispersive,  in  apparent  disagreement  with 
existing  models  of  resonant  nonlinear  response.®  We 
confirm  this  general  feature  of  avalanche  nonlineari¬ 
ties  that  nonresonant  contributions  dominate  the  non¬ 
linear  polarization  of  this  type  of  resonant  interaction 
and  point  out  several  important  applications. 

Basic  features  of  avalanche  absorption  were  mod¬ 
eled  by  Kueny  et  alJ  and  others.®  Those  authors 
pictured  the  cross  relaxation  of  one  excited-  and  one 
ground-state  ion,  which  resulted  in  two  excited- 
intermediate-state  ions,  as  a  key  step  in  the  under¬ 
lying  nonlinear  process  leading  to  a  chain  reaction  of 
excitation  driven  by  excited-state  absorption.  Ni®  sub¬ 
sequently  pointed  out  the  essential  role  of  migration 
in  the  process,  and  numerous  examples  of  avalanche 
absorption  have  since  been  reported.^®  Significant 
applications  have  emerged  in  the  areas  of  avalanche 
upconversion  lasers,^  amplifiers, and  signal  process¬ 
ing.^®  However,  there  have  been  few  fundamental 
studies  of  the  nonlinear-optical  response  associated 
with  the  avalanche  itself.  Here  we  provide  direct 
evidence  that  the  real  part  of  the  nonlinear  index 
is  much  larger  than  the  imaginary  part,  exactly  in 
resonance  with  the  avalanche  transition. 

In  conventional  theories  of  nonlinear  optics  a  single 
strong  input  field  in  resonance  with  an  allowed  tran¬ 
sition  generates  a  purely  absorptive  change  in  the  re¬ 
sponse.  On  the  other  hand,  at  large  detunings  the 
response  is  predicted  to  be  dispersive.  In  avalanche 
upconversion  an  interesting  situation  is  encountered  in 
which  two  transitions  are  involved,  one  exactly  on  reso¬ 
nance  and  one  far  off  resonance,  as  indicated  at  the  left 
in  Fig,  1(a).  Intuitively,  in  the  absence  of  coupling  be¬ 
tween  the  atoms  one  expects  the  resonant  interaction 
to  dominate.  Hence  the  nonlinear  response  should  be 
absorptive  in  character. 

For  coupled  atoms,  though,  the  situation  is  greatly 
modified.  Above  avalanche  threshold,  cross  relax¬ 
ation  provides  a  mechanism  for  heavily  populating 


onversion  transition  in  concentrated  Tm  LAT;  confirm 
i  optical  interaction,  the  induced  response  is  strongly 
be  a  general  feature  of  avalanche  polarization,  with  an 
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excited  states  even  when  the  initial  optical  interaction 
involving  the  ground  state  is  far  off  resonance.  The 
populations  in  excited  states  and  the  value  of  the  pre¬ 
dicted  threshold  intensity  show  little  or  no  dependence 
on  the  weak  ground-state  absorption.'^  However, 
models  that  omit  the  ground-state  transition  alto¬ 
gether  predict  only  nonlinear  properties  associated 
with  the  excited-state  transition,  which  is  purely 
absorptive  on  resonance.  Because  the  excited  state 
involved  in  the  off-resonant  ground-state  transition 
can  nevertheless  be  significantly  populated  because  of 
avalanche  processes,  the  polarization  associated  with 
the  ground-state  optical  interaction  is  unexpectedly 
strong  even  when  the  ground-state  absorption  is 
negligible. 


Fig.  1.  (a)  Reduction  of  a  multilevel  avalanche  system  to 
a  two-level  nonlinear  system.  Cross  relaxation  (curved 
arrows)  and  decay  (wiggly  arrows)  account  for  incoher¬ 
ent  pumping  of  excited  states  at  large  detunings  from  the 
ground-state  transition,  (b)  Calculation  of  the  suscepti¬ 
bility  from  Eq.  (4).  Solid  dashed  curves  correspond  to  real 
and  imaginary  induced  indices,  respectively.  The  parame¬ 
ters  used  are  r2i  ==  2  x  10^^  s“^  and  721  =  2  x  10® 
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To  model  this  situation  we  replaced  all  but  one  ex¬ 
cited  state  with  an  avalanche  reservoir  whose  occupa- 
tion  IS  large  but  bears  no  particular  phase  relationship 
to  the  polarization  developed  on  the  off-resonant  tran¬ 
sition  between  the  ground  state  and  level  12).  This 
quasi-two-level  model  is  shown  at  the  right  in  Fig.  lia). 
Incoherent  pumping  and  depletion  of  the  reservoir  are 
justified  in  the  limit  of  large  A  by  the  spontaneous  na¬ 
ture  of  the  cross  relaxation  that  populates  state  |2>  and 
by  spontaneous  decay  processes  among  excited  levels. 
In  steady  state  the  equations  of  motion  for  the  compo¬ 
nents  of  the  density  matrix  are 

Pii  =  y2iP22  -  Ht’ - -\i  =  0,  (1) 

P22  ~  ~y-l\PT2,  fit’  .\2  =  0  ,  (2) 

P21  “  ~ir2i  -  iA)p2i  P22  —  Pll)  =  0.  (3) 

Here  ■>'21  and  r2i  are  the  population  decay  and  the 
optical  dephasing  rates,  respectively.  A  =  w  -  W21 
is  the  detuning  of  the  incident  frequency  <v  from 
transition  frequency  C021  between  levels  |1)  and  |2),  H  = 
-2fX32E/h  for  fields  of  the  form  E  =  E  exp(-iwi)  -t- 
C.C.,  P21  is  the  slowly  var5dng  amplitude  of  the  off- 
diagonal  element  P2i  defined  by  p2i  =  P21  exp(-:a)<), 
and  V  is  the  imaginary  part  of  pi2  =  P21  =  u  + 

IV.  Closure  requires  that  pu  -t-  P22  +  pr  =  1,  where 
reservoir  population  pr  gives  the  total  population 
density  circulating  among  excited  states  other  than  |2>. 
In  steady  state  .\i  =  A2,  and  we  assume  that  this  rate 
is  proportional  to  occupation  of  the  reservoir  (A2  =  f  pr, 
where  f  is  the  rate  at  which  reservoir  atoms  decay  to 
level  |2)). 

The  steady-state  polarization  P(t)  =  N  Tripp.)  = 
.P(w)exp(  +  c.c.  can  be  determined  from  Eqs.  (D- 
(3)  and  used  to  calculate  an  effective  susceptibility 
^eff  according  to  P{aj)  —  Xefr  consists  of 

the  usual  linear  susceptibility  and  a  nonlinear, 

intensity-dependent  part  xl^ico).  The  nonlinear  part 
is^"* 


-  4iV|pi2l^/i~^€o-i(A  -  ir2i)r2i/y2i 
.(A2  +  r2i2)(A2  r2i2  +  mPTai/yzi) 


4iV|pi2pA~No-HA  -  ir2i)(l  +  2f/y2i) 

A2  +  r2i2  +  ini2r2i/y2i 


(Pr/|fl|2)  . 


(4) 

The  first  term  in  Eq.  (4)  corresponds  to  saturation 
of  a  conventional  two-level  system,  whereas  the  sec¬ 
ond  term  arises  from  the  incoherent  pumping  channel. 
The  relative  contributions  of  the  two  terms  for  vari¬ 
ous  detunings  of  the  light  wave  from 

the  excited- state  transition  frequency  vary-  consider¬ 
ably.  For  small  detunings  [A'  =  d  -  (wba  -  ca2i)  ^  0] 
the  first  term  may  be  much  smaller  than  the  sec¬ 
ond,  because  detuning  from  the  ground-state  transi¬ 
tion  under  this  condition  (A  =  cjfja  ~  (^21)  is  necessarily 
large  for  avalanche  processes.  In  this  case  the  nonlin¬ 
ear  susceptibility  is  dominated  by  a  real  contribution 
that  peaks  at  the  excited-state  resonance,  as  shown 
in  Fig.  Kb)  (A'  =  0  ^  A  =  20).  On  the  other  hand, 
when  A  =  0  the  first  term  can  dominate,  because  pr 


tends  to  zero  abruptly  lA  0)  and  the  avalanche  is 
extinguished.  This  situation  generates  broad  features 
in  the  vicinity  of  the  ground-state  resonance  [A  =  0 
in  Fig.  I'b)],  where  the  real  nonlinear  index  passes 
through  zero. 

The  total  susceptibility  is  the  sum  of  the  off-resonant 
ground-state  polarization  from  Eq.  (4)  and  the  reso¬ 
nant  excited-state  polarization  that  one  obtains  by 
ignoring  the  ground-state  transition.  The  anaUtical 
results  in  Fig.  Kb)  for  optical  susceptibility 
based  on  microscopic  avalanche  dynamics  can  be  com¬ 
pared  directly  with  measurable  index  changes  ob¬ 
served  by  two-beam  coupling.  In  the  weak  probe  limit, 
experimental  values  for  the  Debye  decay  time  r  that 
characterize  the  saturation  process,  as  well  as  real  and 
imaginary  parts  112'  and  rto"  of  the  nonlinear  refractive 
index  1x2  (directly  proportional  to  4eff  iK  can  be  deter¬ 
mined  from  the  analysis  of  TBC  spectra. 

Experiments  were  performed  with  a  pump-probe 
technique  that  we  described  previously,  which  makes 
use  of  a  stabilized  dye  laser.^  Only  one  improvement 
over  our  earlier  approach  was  implemented  here.  The 
output  was  split  into  two  beams  propagating  with  po¬ 
larizations  parallel  to  the  optic  axis  and  a  third  beam, 
which  did  not  pass  through  the  sample  but  provided  an 
intensity  reference  for  a  differencing  signal  amplifier 


-0.3H 


Fig.  2.  Experimental  two-beam  coupling  signals  in 
Tm:LiYF4  versus  both  pump -probe  detuning  and  pump 
intensity. 


Fig.  3.  Real  (^2')  2ind  imaginary  (n2'0  nonlinear  indices 
versus  input  intensity  from  data  of  Fig.  2. 
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Wavelength  (nm) 

Fig.  4.  Plot  of  the  dispersive  (filled  diamonds)  and  ab- 
amplitudes  of  the  TBC  signals  in 
o/c  lm.LiYf4  versus  wavelength  over  a  broad  spectral  re¬ 
gion  containing  both  ground-  and  excited-state  transitions. 
Inset.  Enlargement  of  the  avalanche  resonance  region  at 
8  nm.  The  slight  skewing  of  the  dispersive  component 
line  shape  with  respect  to  the  absorptive  component  is 
due  to  the  summation  of  ground-  and  excited-state  po- 
al)s*(D^t°*^  Data  are  uncorrected  for  sample 


to  reduce  common  mode  intensity  noise  by  40  dB.  The 
samples  were  single  crystals  of  Tm;LiYF4  with  concen¬ 
trations  of  1  - 12%.  The  strongest  avalanche  response 
tree  of  thermal  grating  contributions  was  obtained  in 
a  1.5%  sample  5  mm  in  length,  which  we  used  to  mea¬ 
sure  the  intensity  dependence  reported  here.  A  5% 
sample  3  mm  in  length  was  used  for  all  the  wide- 
spectral-range  measurements.  Typical  spectra  re¬ 
corded  at  various  intensities  are  shown  in  Fig  2  A 
strildng  feature  of  these  data  is  the  rapid  rise  in 
nonlinear  response  above  a  well-defined  threshold 
intensity  near  20  kW/cm^.  The  shapes  of  individual 
spectral  traces  are  well  described  by  the  superposition 
ot  absorptive  and  dispersive  components  for  a  single 
nonlinear  resonance.^ 

The  values  of  112'  and  na"  determined  from  TBC 
analysis  are  shown  in  Fig.  3.  An  abrupt  rise  in  nonlin¬ 
ear  response  occurs  just  above  the  avalanche  threshold 
at  20  kW/cm  .  Then,  in  the  range  40-100  kW/cm^, 
both  the  real  and  the  imaginary  parts  diminish  slowly 
because  of  saturation.*^  The  most  interesting  aspect 
of  these  data  is  that  the  real  part  m'  of  the  induced  in¬ 
dex  is  much  larger  than  the  imaginary  part  n2".  This 
behavior  is  contrary  to  predictions  of  resonant  response 
in  a  two-level  system. 

These  characteristics  are  nevertheless  exactly  these 
predicted  by  Eq.  (4)  when  processes  such  as  cross 
relaxation  exist  that  permit  incoherent  pumping  of  an 
excited  state  in  the  presence  of  off-resonant  excitation 
of  ground-state  centers.  Figure  1(a)  clearly  shows 
the  change  in  nonlinear  response  predicted  by  our 
model  when  incident  light  is  tuned  to  a  ground-state 
absorption  resonance  (A  =  0)  or  to  an  excited-state 
resonance  (A'  =  0  A  =  20).  These  results  are  in 
excellent  agreement  with  measurements  (Fig.  4)  that 
show  that  712'  is  zero  at  the  686.5  nm  ground-state 
resonance  but  nonzero  and  dominant  at  the  648.2-nm 
excited-state  (avalanche)  resonance. 


Our  result.s  show  that  the  avalanche  nonlineantv 
obser%-ed  in  Tm:Li\T4  at  648.2  nm  gives  rise  to  beam 
coupling  with  a  sharp  threshold  and  a  dispersive 
character  at  room  temperature.  The  real  and  the 
imaginary  parts  of  the  induced  nonlinear  index  have 
large  magnitudes,  of  the  order  of  10  **  m^  \V.  The 
real  index  pes  to  zero  as  expected  at  ground-state  reso¬ 
nant  wavelengths  but  is  large  and  dominant  over  the 
imaginaiyv  index  at  the  excited-state  ( avalanche  (  reso¬ 
nance.  This  behavior  is  consistent  with  a  simple 
density-matrix  model,  which  reveals  that  the  dis¬ 
persive  character  of  the  nonlinear  susceptibility  is 
a  general  feature  of  any  system  that  experiences 
significant  pumping  of  an  excited  state  at  large 
detunings  from  any  ground-state  resonance.  This 
dispersive  nonlinearity  could  provide  a  new  passive 
self-mode-locking  mechanism  for  avalanche  upconver- 
sion  lasers  and  presents  an  intrinsic  limitation  to  gain 
in  avalanche  upconversion  fiber  amplifiers. 

H.  Ni  thanks  the  Barbour  Scholarship  Program 
for  graduate  support.  The  authors  are  indebted  to 
N.  Cockcroft  and  R.  M.  MacFarlane  for  samples. 

References 

1.  M.  A.  Kramer,  W.  R.  Tompkin,  and  R.  W.  Boyd,  Phvs. 

^  <1986);  I.  McMichael,  P.  Yeh,  and 

P.  Beckwith,  Opt.  Lett.  13, 500  ( 1988). 

2.  S.  A.  Boothroyd,  J.  Chrostowski,  and  M.  S.  O’Sullivan 

J.  Opt.  Soc.  Am.  B  6,  766  (1989). 

?7,^n5"(  1992) 

124  a97T)’ 

t’  “Aval^che  upconversion  in  Tm;LiYF4, 

TmiYAG,”  Ph.D.  dissertation  (Univer¬ 
sity  of  Michigan,  Ann  Arbor,  Mich.,  1994) 

6.  See,  for  example,  Y.  R.  Shen,  The  Principles  ofNonlin- 
ear  Optics  (Wiley,  New  York,  1984) 

’  B  ft  SS  ^  ^ 

A.  W.  Kueny,  J.  Lumin. 
M.  J.  Riley,  P.  S.  May, 
and  H.  U.  Gudel,  Coord.  Chem.  Rev.  Ill,  125  (1991)’ 

« lorafSai  ““  ^  ® 

io  1424  (1991). 

10.  T  Hebert,  R.  Wannemacher,  R.  M.  Macfarlane.  and 
W  Lenth,  Appl.  Phys.  Lett.  60,  2592  (1992);  Y.  H. 
Chen  and  F.  Auzel,  Electron.  Lett.  30,  1602  (1994)-  M 
Pollnau,  W.  Luthy,  and  H.  P.  Weber,  J.  Appl.  Phys  V?' 
®J28  (1995);  R.  Scheps,  IEEE  J.  Quantum  Electron.  31, 

11.  T.  San^ock,  E.  Heumann,  and  G.  Huber,  in  Interna- 
tional  Quantum  Electronics  Conference  (Optical  Society 
of  America,  Washington,  D.C.,  1996). 

12.  L.  M.  Yang,  T.  Sonsnowski,  M.  L.  Stock,  T.  B.  Norris 

J.  Squier,  G.  Mourou,  M.  L.  Dennis,  and  I.  Duling,  Ont’ 
Lett.  20, 1044  (1995).  ^ 

13.  H.  Ni  and  p.  C.  Rand,  Opt.  Lett.  17,  1222  (1992). 

14.  Q.  Shu,  Cooperative  optical  nonlinearities,”  Ph.D. 
dissertation  (University  of  Michigan,  Ann  Arbor,  Mich.^ 

15.  pr  can  be  calculated  by  use  of  the  standard  four-level 
model,  as  we  have  done,  or  it  can  be  treated  as  a  free 
parameter. 


PHYSICAL  REVIEW  B 


VOLUME  55.  NUMBER  14 


1  APRIL  1997-n 


Critical  slowing  down  and  dispersion  of  avalanche  upconversion  dynamics 

Q.  Shu  and  S.  C.  Rand 

Division  of  Applied  Physics.  1049  Randall  Laboratory.  University  of  Michigan.  Ann  Arbor.  Michigan  48109-1120 

(Received  23  September  1996) 

The  temporal  response  time  on  an  avalanche  upconversion  transition  is  shown  to  undergo  critical  slowing 
down  near  the  avalanche  threshold.  Also,  the  nonlinear  refractive  index  is  predominantly  real,  despite  the 
resonant  nature  of  the  avalanche  phenomenon.  These  findings  are  in  excellent  agreement  with  density-mauix 
analyses  showing  that  cross  relaxation  dramatically  alters  the  effective  excited-state  lifetime  and  that  off- 
resonance  absorption  curiously  dominates  the  wavelength  dependence  of  the  nonlinear  refractive  index. 
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I.  INTRODUCTION 

Since  its  discovery,'  avalanche  upconversion  has  been  of 
interest  for  solid-state  lasers^  because  of  the  strong  absorp¬ 
tion  it  generates  in  spectral  regions  normally  devoid  of 
ground-state  optical  transitions  in  concentrated  rare-earth 
and  transition-metal  crystals.  Materials  displaying  the  ava¬ 
lanche  phenomenon  universally  exhibit  cross-relaxation  dy¬ 
namics  in  connection  with  an  excited-state  optical  resonance. 
By  tuning  light  into  coincidence  with  an  excited-state  ab¬ 
sorption  resonance  and  providing  sufficient  power  density  to 
exceed  the  avalanche  threshold,  intense  luminescence  can  be 
generated  at  wavelengths  shorter  than  the  incident  one. 

Most  earlier  investigations  focused  on  population  inver¬ 
sion  conditions  in  avalanche  systems.  Only  a  few  studies 
have  examined  the  nonlinear  response^'  associated  with  ava¬ 
lanche  dynamics.  In  an  earlier  publication,*  we  reported 
measurements  of  nonlinear  refractive  indices  as  a  function  of 
incident  intensity,  using  a  two-beam  coupling  approach. 
Here  we  describe  experimental  observations  of  a  previously 
unreponed  characteristic*  of  avalanche  nonlinearities, 
namely  critical  slowing  down  of  the  effective  response  time. 
We  furnish  the  theoretical  basis  necessary  to  understand  this 
result,  together  with  a  detailed  explanation  of  the  frequency 
dependence  of  the  nonlinear  refractive  index.  Avalanche  up¬ 
conversion  is  shown  to  be  an  example  of  surprising,  but  quite 
general,  dynamics  in  which  the  refractive  index  changes  as¬ 
sociated  with  polarization  on  an  excited-state  resonance  are 
primarily  real,  controlled  by  the  off-resonant  contributions 
from  weak,  ground-state  absorption. 

n.  THEORY 
A.  Steady-sUte  dynamics 

To  formulate  time-dependent  analysis  of  dynamics  in  an 
avalanche  system,  it  is  first  necessary  to  establish  a  steady- 
state  description  of  the  occupation  of  energy  levels.  To  this 
end  we  begin  by  considering  the  atomistic  model  of  die  ava¬ 
lanche  process  illustrated  in  Fig.  1.  The  basic  dynamics  con¬ 
sist  of  excitation  on  an  excited-state  transition,  followed  by 
internal  relaxation  to  a  level  which  participates  in  cross  re¬ 
laxation  with  a  second  atom  in  its  ground  state.  The  Hamil¬ 
tonian  for  such  a  system  can  be  written  as  H  — 
where  Hq  is  the  unperturbed  Hamiltonian  and  describes 
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its  interaction  with  the  field  (W,iom-fieid='“M  E)  and  with 
other  atoms  Treating  relaxation  processes  other 

than  cross  relaxation  phenomenologically,  the  equation  of 
motion  for  the  density  matrix  is 

1  ftp =[W,p]-i- relaxation  terms.  (1) 

Hence  equations  for  the  temporal  evolution  of  individual  el¬ 


ements  of  the  density  matrix  are 

Pll=  r2P22+  r3lP33+  y*lPu- aPuP33 

-  r2P22+  y^2P^i'^  y42P44  +  2apiiP33 

-2^p|2-iV24P42+jV42P24.  (3) 

P33=  -  r3P33+  r43P44-  “Pi  lP33+  ^pL'*’  ^^^Pl  I  “  P33). 

(4) 

P44=-y4P44+tV'42P24"‘'^24P42-  (3) 

P24~('“42“r42)P24~  <^24(P44“P22)> 


FIG.  1.  Schematic  illustration  of  the  basic  dynamics  in  the  ava¬ 
lanche  model.  Absorption  of  light  on  the  transition  between  states 
|2)  and  \4)  is  followed  by  internal  relaxation  to  state  |3)  from  which 
cross  relaxation  with  a  ground-state  neighbor  can  occur. 
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P42“P24'  (7) 

In  these  equations,  0)42  is  the  optical  transition  frequency 
t>etween  state  |.2)  and  l4>,  the  y,  is  the  population  decay  rate 
of  level  i  (i  =  2,3,4),  r42  is  the  dephasing  rate  of  the  optical 
transition  (including  dephasing  caused  by  energy  migration), 
a  is  the  cross-relaxation  coefficient,  and  /3  is  the  coefficient 
for  the  reverse  reaction.  B  is  the  Einstein  coefficient  of 
stimulated  absorption  and  emission,  and  I  is  the  input  inten¬ 
sity* 

Taking  the  input  field  to  be  of  the  form  E-Ee  '‘^'4-c.c., 
the  matrix  element  of  the  atom-field  interaction  in  the  rotat¬ 
ing  wave  approximation  can  be  written  as 

^24=  V^4*2=(2|//aton..fieldl4)=  ^  024^'"'.  (8) 

The  quantity  0.24= -2ii24E*/h  is  related  to  the  intensity  ac¬ 
cording  to 


sity  changes  very  little.  In  the  limit  g=0,  the  occupation  of 
the  excited  state  can  be  written  as 


Clearly,  for  aj<0  one^  obtains  P22=0.  whereas  for  ai>0 
one  finds  the  nonzero  value  P22““^i/<32’  the  excited 
state  becomes  occupied.  In  this  limit  the  boundary  of  the 
dynamical  phase  transition  which  occurs  at  avalanche  thresh¬ 
old  is  well  defined  and  given  simply  by  the  condition  ai=0. 
The  threshold  does  not  depend  sensitively  on  detuning  from 
the  ground-state  transition  under  these  conditions,  but  does 
depend  strongly  on  detuning  from  the  excited-state  reso¬ 
nance.  When  dephasing  and  detuning  are  included,  the 
threshold  condition  for  avalanche  occurrence  can  be  given  in 
terms  of  a  threshold  value  for  /  which  is 


1 


ech 


/=-ect2£|2  =  ^|n24i 

4  ^^24 


(9) 


_ riia+y-j) _ 

(“-  r3i)(r43-  r4i)-(r3  +  r3i)r4i ' 


(20) 


where  iJi^i={2\-et\A)  is  the  transition  dipole  moment,  e  is 
the  dielectric  constant  of  the  medium,  c  is  the  speed  of  light 
in  vacuum,  and  h  is  Planck’s  constant. 

In  a  closed  system  we  can  easily  solve  for  steady-state 
behavior.  The  results  are 


P44==/p22> 

(10) 

P33“^P22“*“ 

(11) 

-a^-yja]- 40200 

(12) 

202 

The  ground-state  occupation  is  determined  by  the  closure 
relation  Pn+P22'*‘P33'^P44~l*  For  convenience  in  simplify¬ 
ing  these  expressions,  several  parameters  have  been  intro¬ 
duced: 


/- 


l/2|n24|^r42/ ^4 


k=-g 


A2-fn2+i/2in24|T42/r4’ 

4fl'|n24l^ 

«  r3+y3,  +  4fi'lft24|^’ 

1+/  ,  (y43-r4i)/-r2 


2  73  +  731 

ao=g[a(l-  i^)+  J  732]. 


(1-g). 


(13) 

(14) 

(15) 

(16) 


ai=  -72~(74i+  743)/+[732+2a(l-g)]*-ag(l+/  ), 

(17) 

02= -2a(l +/+<:)* -2y8,  (18) 

The  quantity  5'  in  Eq.  (14)  is  a  scaled  Einstein  coefficient 
given  by  B'  =  Bll\024\'^ . 

Notice  that  the  quantity  g  is  related  to  the  off-resonance 
absorption  rate.  As  g  approaches  zero,  the  appearance  of 
avalanche  emission  above  a  certain  intensity  becomes  more 
and  more  abrupt,  although  the  value  of  the  threshold  inten- 


To  complete  these  results  we  add  the  solution  for  the  off- 
diagonal  element  P24,  which  is  of  the  form  P24=^24^''^^  P24 
being  its  slowly  varying  amplitude.  Writing  out  this  ampli¬ 
tude  explicitly  in  terms  of  real  and  imaginary  parts  according 
to  P24=u-^iv,  the  expectation  value  of  the  atomic  polariza¬ 
tion  can  then  be  determined  from  the  relation 

p(r)  =  Tr(p/r)  =  /X24(w-/i;)e“'^'  +  c.c.,  (21) 

for  wavelengths  in  the  vicinity  of  the  avalanche  transition. 
The  solutions  for  u  and  v  are  given  by 


474  , 

~  n  r  -^^22  ’ 

4*24*  42 

(22) 

y  =  7r-/P22- 
4*24 

(23) 

The  macroscopic  nonlinear  polarization  density  Pit) 
=  Np(0  has  a  slowly  varying  amplitude  given  by 
P{o))-N fi2A^u-iv)  from  which  the  effective  nonlinear 
susceptibility  can  be  determined.  This  is  accomplished  by 
first  identifying  the^ull  nonlinear  susc^tibility  accord¬ 
ing  to  the  relation  P^^'(a>)=;^^^^^(cu)E(c*>),  and  then  picking 
out  the  effective  susceptibility  using  the  correspon¬ 

dence  ;^^^^(cu)E(a>)=;^eff|£(ct>)l%(cu).  The  resulting  ex¬ 
pression  for  the  effective  third-order  susceptibility  is 


2T2,\E\* 


(4-1T24). 


(24) 


Notice  that  the  real  part  of  this  expression  is  zero  exactly 
on  resonance  (A=0).  This  is  the  usual  expectation  for  optical 
response  in  two-level  systems.^ 


B.  Time-dependent  dynamics 

For  an  avalanche  system  which  is  not  in  a  steady  state,  the 
nonlinear  nature  of  the  dynamical  equations  suggests  that 
numerical  techniques  must  be  used  to  solve  for  system  evo- 
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RG.  2,  Phase  diagram  of  temporal  evolution  in  a  nonlinear  dy¬ 
namic  system,  illustrating  the  idea  of  linearization  of  the  equation 
of  motion  for  the  density.  matri;c  p.  Initially  the  system  resides  in 
state  /4,  but  is  perturbed  by  an  intensity  change  which  causes  to 
become  nonzero.  The  system  subsequently  evolves  to  state  C,  for 
which  point  D  is  an  adequate  approximation  for  small  A/. 

lution.  However  when  changes  are  slow  or  small,  an  analytic 
approach  is  still  possible,  based  on  linearization  of  the  dy¬ 
namical  equations.^ 

The  basis  for  the  linearization  method  is  illustrated  in  Fig. 
2.  This  diagram  shows  a  phase  space  of  solutions  to 
the  equations  of  motion  of  the  density  matrix  for  two  input 
intensities  /  and  /-+- A.  Initially,  the  system  resides  in  steady 
state  A.  When  a  small  intensity  step  A/  is  applied  to  the 
system  at  f=0,  its  state  is  displaced  to  B  and  subsequently 
follows  a  nonlinear  trajectory  to  reach  a  new  steady  state  C, 
as  time  progresses.  If  A/  is  small  enough,  we  can  approxi¬ 
mate  this  evolution  using  the  tangent  at  B,  Then,  the  line  BD 
replaces-  the  nonlinear  trajectory,  permitting  us  to  find  a  lin¬ 
ear  approximation  for  the  new  steady  state  of  the  system. 
Final  results  were  obtained  by  successive  approximations  in 
the  limit  of  A/--^0. 

In  the  rate  equation  limit,  the  earlier  dynamical  equations 
become 

P22=  “  y2p22+  y32P33+  r42P44+ 2 QfpnPjj- 2/3p22 

“524^(P22“P44)>  (25) 

P33=  “  r3P33+  r43P44-“PnP33  +  /3p^  +  fl|3/(Pll-p33). 

(26) 

P44=-y4P44+fl24/(P22-P44).  (27) 

5  24  and  fl  13  are  the  effective  pumping  rates  on  transitions 

\2)—^\4)  and  |l>“^|3),  respectively. 

Locally  linearized  versions  of  these  equations^  are  ob¬ 
tained  with  the  following  substitutions: 

P„-[p„  +  ^,(0)]/2,  (28) 

P„P;V^[P/.Py7(0)  +  P»(0)P;;]/2,  (29) 

P„-[p„  +  P„(0)]/2  (30) 

Indices  ij  can  both  take  on  the  values  2,  3,  or  4.  After 
substitution,  the  equations  become 

(31) 


P33~^32P22'''^33P33‘''‘^34P44'*"<^3>  (32) 

P44=a42p22+U43P33"'"^44P44'*'<^4-  (33) 

The  coefficients  are  defined  as  follows: 

^22~  ~[y2'*'®247"*"4^P22(0)  +  2ap33(0)],  (34) 

^23=y32  +  2a-2Q;p22(0)-4ap33(0)-2ap44(0), 


(35) 

a24=y42  +  524/-2ap33(0),  (36) 

^2=  “(y2  +  S24/)p22(0)  +  (y32+2ar)p33(0) 

+  (y42+fl24/)P44(0)-p22(0),  (37) 

£2  32=  “  l3^~2/3p22(0)  -  ap33(0)],  (38) 

2233=  -  [  y3  +  2S 13/  +  a-  ap22(0)  -  2ap33(0)-  ap44(0)], 

(39) 

yi3~ apjj{Q),  (40) 

C3  =  25i3/— Si3/p22(0)  — (■y3  +  2fl|3/+a)p33(0) 

+  (y43-S,3/)P44(0)-p33(0),  (41) 

£242=524/.  (42) 

£243  =  0,  (43) 

£244=-y4-524/.  (44) 

C4=524/P22(0)~(y4  +  524/)p44(0)-p44(0).  ,  (45) 


Initial  population  densities  required  for  the  coefficients  in 
Eqs.  (3l)-(33)  are  determined  by  steady-state  values  ac¬ 
quired  prior  to  r =0  as  determined  in  Sec.  II  A,  for  a  specified 
intensity  level.  Initial  time  derivatives  are  calculated  by  sub¬ 
stituting  these  initial  populations  directly  into  Eqs.  (25)- 
(27). 

The  equations  of  motion  can  now  be  written  in  a  conve¬ 
nient  matrix  form: 

p=Ap4-c,  (46) 

Here 

/  P22\  /  ^22  ^23  ^24\  /  ^2\ 

p=|  P33  L  A=  ^32  ^33  ^34  »  C=  ^3!. 

\  P44/  \  ^42  ^43  ^44/  \cil 

Applying  the  method  of  Laplace  transforms,  and  represent¬ 
ing  the  transform  of  p  by  R(5),  the  equation  of  motion  yields 

1 

^R(j)-p(0)=AR(5)+ jc.  (47) 

The  Laplace  transform  of  the  solution  is  therefore 

J  1 

R(i)  =  (S-/l)-‘  p(0)+-c  .  (48) 

s 

where  5  is  5  times  the  unit  matrix  and  p(0)  is  the  initial  value 
of  p. 


P 22  “  ^  22P22  ^  23P33 ^ 24P44  ^  ^  2  » 
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FIG.  3.  The  effective  relaxation  time  of  state  |2)  versus  input 
intensity.  The  population  difference  ip22~Pu)  is  graphed  for  direct 
comparison  of  the  peak  in  the  time  constant  with  the  threshold  of 
avalanche  absorption.  Inset:  Divergence  of  the  response  lime  r-  as 
the  ground-state  absorption  rate  5,3  approaches  zero,  calculated 
numerically. 


R(.v)  may  be  assumed  to  consist  of  a  transient  pan  and  a 
steady-state  pan.  That  is,  R(s)  =  R^\s)~^R''\s).  The  steady- 
state  pan  is  readily  identified  as  R'"(.v)=  -( l/5)/\"'c  since 
p{^)  =  lim,_,)[.i'R(.v)]=-/\“‘c.  The  remainder  of  the  solu¬ 
tion  transform  must  be  the  transient  response.  Consequently, 
the  rih  element  of  this  transient  pan  of  the  vector  R^\s)  is 


.\  = 


R,{s 


m,s~-hn,s  + 1, 


(49) 


where  m, ,  /t, .  /, ,  iVf ,  A/,  and  L  are  all  functions  of  elements 
of  A.  defined  in  the  Appendix. 

In  the  long-time  limit  R^\s)  can  be  analyzed  by 

panial  fractions,  or  by  the  expansion  of  Eq.  (49)  to  first  order 
in  s.  Comparing  the  result 


R'-{s)  = 


LfiiiN/L-nilli) 


(50) 


s-^ 


NIL-n,ll, 


with  the  general  form  of  a  transformed  exponential,  namely 
=  where  c  and  a  are  constants,  we 

find  an  effective  time  constant  of 


ii 


(51) 


for  level  /.  Calculations  of  the  effective  time  constant  Ti  (for 
level  2)  versus  intensity  in  Fig.  3  reveal  a  dramatic  “slowing 
down”  at  the  avalanche  threshold.  The  inset  of  Fig.  3  shows 
that  Ti  in  fact  diverges  at  a  critical  intensity  as  the  ground- 
state  absorption  rate  5 13  approaches  zero.  This  is  a  key  result 
of  the  present  research. 


C.  Wavelength  dependence 

A  significant  limitation  of  conventional  analyses  of  the 
avalanche  effect  is  that  off-resonant  excitation  from  the 
ground  state  which  provides  the  seed  population  in  the  ab¬ 
sorptive  excited  state  is  ignored.  This  omission  is  justified 
when  steady-state  populations  are  the  only  quantities  of  in- 


FIG.  4.  A  simplified  model  for  the  description  of  complex  non¬ 
linear  susceptibility  in  avalanche  systems.  On  the  left  is  a  multilevel 
sketch  refecting  details  of  avalanche  dynamics.  On  the  right  is  a 
quasi-two-level  model  which  retains  aspects  essential  for  under¬ 
standing  the  nonlinear  dispersion  of  avalanche  nonlinearities. 

terest,  for  they  depend  negligibly  on  initial  conditions  (such 
as  excited-state  occupation).  However  to  understand  the 
wavelength  dependence  of  the  nonlinear  refractive  indices, 
the  nonlinear  dispersion,  it  is  necessary  to  account  for  the 
initial  (ground-state)  absorption  step  in  the  avalanche  process 
because  of  a  surprising  enhancement  of  the  off-resonant 
ground-state  polarization  by  the  large  excited-state  popula¬ 
tions  generated  through  avalanche  cross  relaxation. 

In  this  section,  the  important  role  of  the  off-resonant  ab- 
.sorption  from  the  ground  state  which  initiates  the  avalanche 
in  determining  polarization  properties  of  the  system  is  exam¬ 
ined.  Rather  than  calculating  populations  and  polarizations 
from  a  five-level  model  which  incorporates  two  optical  tran¬ 
sitions  explicitly,  it  is  very  helpful  to  solve  for  the  nonlinear 
dispersion  in  an  equivalent  two-step  fashion  using  a  quasi- 
three-level  model  which  renders  the  basic  physics  much 
more  obvious.  To  highlight  the  unexpectedly  important  role 
of  the  ground-state  transition  in  determining  nonlinear  dis¬ 
persion,  we  reduce  the  basic  avalanche  model  to  that  shown 
on  the  right  side  of  Fig.  4,  in  which  all  excited  states  but  one 
are  replaced  by  an  incoherently  pumped  reservoir.  The  inco¬ 
herent  pumping  rate  Aj  is  due  exclusively  to  the  nonradiative 
process  of  cross  relaxation,  whereas  A2  contains  contribu¬ 
tions  from  both  cross  relaxation  and  natural  decay  between 
levels.  Both  are  functions  of  intensity.  We  show  that  detun¬ 
ing  from  the  ground-state  transition  determines  whether  the 
nonlinear  response  is  absorptive  or  dispersive. 

After  simplifying  Eqs,  (2)-(7)  to  reflect  the  model  in  Fig. 
4,  the  equations  of  motion  yield 


Pll=r2lP22-^^U-.\,=0,  ,  (52) 

P22=“r2lP22  +  i^l'''’A2  =  0,  (53) 

P21=~(r21-lA)p21+  5'f^(P22“Pll)  =  0-  (54) 

A  =  a>“a>2i  is  the  detuning  of  incident  light  at  frequency  o) 
from  the  transition  frequency  ct>2i  between  levels.  The  clo¬ 
sure  relation  Pii+P22“^Pr=l  accounts  for  the  population 
density  p^  circulating  in  the  avalanche  reservoir.  The  reser¬ 
voir  population  p^  can  serve  as  a  free  parameter  or  may  be 
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determined  explicitly  from  a  multilevel  model  which  takes 
the  second  photon  interaction  and  cross  relaxation  into  ac¬ 
count.  That  is 


P 


P 


a ' 


where  /  runs  over  the  indices  of  levels  other  than  |1)  and  [2) 
in  the  full,  microscopic  model. 

Using  R  to  denote  the  rate  per  atom  of  decay  to  level  |2) 
from  the  reservoir,  we  can  apply  the  steady-state  relation 
Ai  =\^  =  Rpr  to  solve  Eqs.  (52)-(54)  for  p, , ,  P22,  and  po] : 

i/2ia|-r.|/r2, 

A-^r;,r|n|-r:,/y:, 


( I  R/  y2\)i  ^  l/2|n|‘r2i/y2i 

l/2|0|-r2i/y:i 


Pr^ 


(55) 


(A-+r;|)/?/r:,-  l/2|li|-r,,/y:i 
"  A-+r;,  +  iai-r,,/r:, 

(  1/2)11(A~/T2,)  (l/2)e(A-/r2,) 


(56) 


P:r 


A"+1"2i'^I^^I”^2!^T:i  A"+r2i  +  |il|‘'r2i/y2i 


x|  !  +  —  P,-.  (57) 


The  effective  j^usceptibility  may  be  determined  from 
decomposing  the  expression  for 
into  linear  _and  nonlinear  parts  using 
where  the  linear  portion  is 

given  by  the  usual  expression 

X  (0.)= - . 

the  effective  nonlinear  susceptibility  is  found  to  be 


(A-+ri,){A-+rii+in|'r,,/r.,) 

4N\^^2\*f^-^EoH^-i^'2^)i  1  +2/?/y2i) 

^  FTrTTlWrVrT 


x(p,/|n|') 


(58) 


Notice  that  this  result  is  quite  different  from  the  effective 
susceptibility  in  Eq.  (24),  which  overlooked  the  ground-state 
transition.  At  the  avalanche  resonance,  the  excited-state 
detuning  A'  is  zero.  This  means  A'  =  a>“iu^fl 
=  A-(aj^^-a>2i)  =  0,  with  the  result  that  A  =  a>^^-a>2|.  Un¬ 
der  these  conditions,  the  reservoir  population  is  maximized 
and  the  real  component  of  the  nonlinear  susceptibility  is 
large  on  resonance.  In  fact,  whenever  ^  l^he 

real  pan  exceeds  the  imaginary  part  exactly  on  resonance. 


I 

o 


' 


FIG.  5.  Calculation  of  the  real  and  imaginary  components  of  an 
avalanche-induced  nonlinear  susceptibility.  Parameters  used  in  the 
calculation  were  r2t  =  2Xl{)'-  s"'  and  y2i"2xlO^  s'*  (roughly 
appropriate  for  the  —  transition!,  and  -ciM]  =  2x  !()' ' 
s“‘.  |a|'  =  3xl0‘^ 

A  direct  comparison  of  the  predicted  behavior  of  nonlin¬ 
ear  su.sceptibility  components  at  a  ground-state  re.sonance 
(A=0)  and  at  an  avalanche  transition  (A'  =0)  is  made  in  Fig. 
5.  For  this  plot,  the  reservoir  population  p,.  was  taken  to  have 
the  steady-state  value  determined  from  the  rive-level  model 
in  the  left  side  of  Fig.  4.  Notice  that  the  real  component 
passes  through  zero  for  A=()  (ground-state  transition) 
whereas  it  exhibits  a  maximum  for  A'  =0  (cxcitcd-siatc  reso¬ 
nance).  This  is  a  second  key  result  of  the  present  research. 


III.  EXPERIMENT  AND  RESULTS 

Avalanche  upconversion  can  readily  be  observed  in  1.5% 
Tm:LiYF4  crystals  at  room  temperature  using  dye  laser  ex¬ 
citation  at  648  nm.*^  In  this  work  we  measured  the  nonlinear 
indices  and  characteristic  response  time  of  the  avalanche 
nonlinearity  in  this  material  by  a  two-beam  coupling  tech¬ 
nique.  Details  of  the  experimental  method,  analysis,  and  val¬ 
ues  of  the  induced  indices  have  been  described  elsewhere.  “ 
Here  we  report  on  the  temporal  response  and  nonlinear  dis¬ 
persion  over  a  wide  wavelength  range.  Both  these  aspects  of 
avalanche  dynamics  show  unique  features  which  provide  an 
important  basis  for  comparison  with  theory. 

Figure  6  is  a  plot  of  the  experimental  response  time  mea¬ 
sured  in  1.5%  Tm:LiYF4  as  a  function  of  incident  intensity. 
The  avalanche  threshold  occurs  at  I  W/cm^  The  important 
thing  to  notice  is  that  both  above  and  below  the  threshold 
intensity,  the  medium  response  is  considerably  faster  than  it 
is  at  the  critical  point.  This  indicates  a  dramatic  slowing 
down  of  the  effective  lifetime  of  the  absorptive  excited  state, 
right  at  the  avalanche  threshold. 

By  performing  ultrahigh  resolution  scans  of  the  beam¬ 
coupling  spectrum  at  many  different  wavelengths  in  the  red 
spectral  region,  additional  results  shown  in  Figs.  7  and  8 
were  obtained.  The  traces  shown  in  Fig.  7  cover  the  region 
from  the  ground-state  resonance  near  688  nm  to  the  excited- 
state  resonance  near  648  nm.  If  is  immediately  evident  from 
the  trend  in  their  shapes,  that  the  refractive  index  changes 
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FIG.  6.  Experimental  response  time  of  the  avalanche  medium 
versus  intensity,  obtained  by  two-beam  coupling. 

induced  by  incident  light  vary  from  primarily  absorptive  at 
688  nm  to  primarily  dispersive  at  648  hm. 

Quantitative  analysis  of  the  real  and  imaginary 
components'*  of  these  and  other  two-beam  coupling  spectra 
in  the  range  640-700  nm  yielded  the  results  shown  in  Fig.  8. 
In  this  hgure  it  is  somewhat  easier  to  see  that  the  real  index 
goes  through  zero  at  the  ground-state  absorption  resonance 
near  688  nm.  By  contrast,  when  the  wavelength  is  tuned  to 
the  avalanche  transition  at  648  nm,  the  real  index  acquires  a 
large  value,  exceeding  the  imaginary  index  by  a  factor  of  2. 


FIG.  7.  Beam  coupling  spectra  at  different  wavelengths  between 
the  ground-state  absorption  at  688  nm  and  the  avalanche  resonance 
at  648  nm. 


FIG.  8.  Decomposition  of  the  nonlinear  susceptibilities  deter¬ 
mined  by  two  beam  coupling  into  real  (diamonds)  and  imaginary 
parts  (open  circles).  The  spectral  range  includes  an  avalanche  reso¬ 
nance  at  648  nm  and  ground-state  absorption  at  688  nm.  The  lower 
trace  gives  the  polarized  absorption  spectrum  {EWC)  taken  at  room 
temperature  which  was  used  to  correct  the  two-beam  coupling  sig¬ 
nal  in  obtaining  the  nonlinear  indices. 


(a) 


Predicted 


(b)  Frequency  Wavelength 


FIG.  9,  (a)  Decomposition  of  the  real  nonlinear  index  (filled 
circles)  into  symmetric  and  antisymmetric  parts.  The  imaginary  in¬ 
dex  is  also  shown  (open  circles),  but  is  well  described  by  a  single 
Gaussian.  Note  the  small  relative  shift  between  the  experimental 
peaks  of  the  two  indices,  (b)  Schematic  diagram  showing  the  quali¬ 
tative  result  of  combining  the  nonlinear  response  associated  with 
both  optical  transitions  giving  rise  to  avalanche  absorption. 
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IV.  DISCUSSION 

The  significant  increase  that  is  observed  in  optical  re¬ 
sponse  time  at  avalanche  threshold  typifies  behavior  ot'  sys¬ 
tems  undergoing  magnetic or  structural phase  transitions. 
It  is  well  known,  for  example,  that  the  characteristic  time  of 
spin  correlations  in  interacting  spin  systems  undergoes  criti¬ 
cal  slowing  down^"^  at  order-disorder  transitions.  The  phase 
transition  in  our  system  consists  of  the  sudden  appearance  of 
an  excited-state  population  above  a  critical  intensity  and  oc¬ 
curs  at  the  threshold  intensity  determined  by  Eq.  (20).  The 
temporal  response  time  does  indeed  diverge  (r— as  the 
ground-state  absorption  rate  approaches  zero.  Hence  we  refer 
to  this  as  critical  slowing  down.  A  spin  analogy  may  be 
particularly  apt,  since  optical  cross  relaxation,  in  which  a 
first  atom  undergoes  nonradiative  deexcitation  while  a  sec¬ 
ond,  coupled  atom  conserves  energy  by  making  a  dipole 
transition  to  an  excited  state,  is  formally  equivalent  to  secu¬ 
lar  spin-spin  interactions.  The  critical  role  of  temperature  in  a 
spin-ordering  transition  is  replaced  here  by  a  critical  depen¬ 
dence  on  optical  intensity.  Certainly  it  would  be  interesting 
to  explore  adaptations  of  the  continuous-spin  Ising  model  to 
avalanche  dynamics.  Also,  additional  experiments  of  the 
type  we  have  reported  here  should  be  very  fruitful  for  more 
detailed  studies  of  avalanche  critical  exponents  and  phenom¬ 
ena. 

Figure  9(a)  shows  the  result  of  a  decomposition  of  the 
measured  indices  into  polarization  contributions  arising  from 


the  two  optical  transitions  in  the  problem.  Figure  9(b)  shows 
schematically  how  the  two  contributions  add  to  give  a 
slightly  asymmetric  real  component  and  an  imaginary  com¬ 
ponent  which  is  diminished  by  an  excited-state  contribution 
of  opposite  sign.  Notice  that  the  decomposition  accounts  sat¬ 
isfactorily  for  the  slight  asymmetry  observed  in  the  real  in¬ 
dex  peak.  The  real  index  peak  is  shifted  with  respect  to  the 
imaginary  index  curve  due  to  a  small  dispersively  shaped 
contribution  from  the  resonant  excited-state  transition. 

These  findings  are  in  excellent  agreement  with  the  quasi- 
two-level  theory  of  Sec.  II  C.  Because  details  of  the  cross¬ 
relaxation  dynamics  are  omitted  in  this  model  of  nonlinear 
dispersion,  it  can  be  concluded  that  the  predominantly  real 
nonlinear  susceptibility  encountered  in  our  experiments  at 
the  excited-state  absorption  resonance  is  not  limited  to  ava¬ 
lanche  systems.  Rather,  it  is  a  propeay  of  all  systems  in 
which  excitation  greatly  detuned  from  any  ground-state  tran¬ 
sition  nevertheless  generates  large  excited-state  populations. 
Avalanche  systems  provide  merely  one  example  where  this 
occurs,  with  the  cross-relaxation  dynamics  constituting  the 
mechanism  whereby  excited  states  develop  large  populations 
nonradiatively. 
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APPENDIX 

Explicit  expressions  for  the  quantities  in  Eq.  (49)  are 

fn2~X  '  [  6/33(^446' 2  "1“  ^ 23^44^' 3  "i"  24^^  33  ^  23^^  34)^  4]  ^  PZ2^  ^ )  ’  (  A  1  ) 

m3  r  ~  '[(a32^244— a42a34)C2  +  24”  ^^22^ 44)^3"^  (^34^22”  ^24^32)^4]  ■•■P33(0).  (A2) 

^4  =  X  ”  ^  [^42^ 33C2  —  « 23^42<^‘3  32^  23  ”  22^^33)^‘4-]  P44{ 0 )  ’  ^  ) 

n2=x"  ‘[(a33a44+ a33fl^-f‘ ^32^23^^44“  ^42^23^2 34+ <^42^24^33 >^■'2”^  23^44(^22^  ^33 ~^<^^44)^^3 

+  ((323^34”  ^24^33)(  ^22“*"^  33 ^44)^4“  (^33’^  ^44)'’^P22(0)]  ■^^^23P33(0)  +  a24P44(^)’  (A4) 

n^  =  X  *  ^32^44)  22'^  ^  ^44)^2  ( ^32^23^44”  ^ 42^  23^  34  ”  ^ 24^ 42^  22 ”  ^ 44^ 42^ 24 ^^22^44"^  ^^44^22)^3 

.  (<324^32”  ^34^22) (^22 ^33“^  <244)t^4]  '^<232P22(0)“  (a22+^44)p33(0)'^<^  34P 44(0)’ 


n4  =  X  *[  — a42a33(^22'^^33"*"  ^44)^2*^  ^23^42(^22"^^33'^  ^44)<^3“*"  (^22^33“^^42^24^33~*”^22^L  <^32^ 23^2  22”  ^33^32^23 

—  23^  34)04]  +  <342P22(^)  ”  (<^22"^  ^33)p44(Q)»  ( A6) 

/2  =  X"  '[(  -a33a^+a44a42a23a34  +  a33a42^2  23^34“^L^42^24”^il^32^23)<^2“^<^23(  “^42^23^34”^44<342^24+^42^24^33^^44^22 
+  ^33^2^)03  +  (^24^33^22“*"  ^23^34^32“*"  ^24^32^23^44”  ^23^34^22^33”  ^23^34^22^44”  ^23^ 34^ 33^ 44 “I"  ^24^33^44 
”^24^33^32^23)^4]  ^^33^ 44P22(^)  ”  ^23^44P33(^)  ^  (^23^34”  ^24^33)  P44(0)’ 
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=  X~  l^cl  22  32^  33  ~  32^^  44^  42^  24  32^  42^  24^^  33  ”  42^^  22^^  33  ~  34^^  42^^  22^  44  “  34^  42^^  33^  44  34^^  42^^  24)  2 

+  (  -  a  -.4^/  42<^  22^  23  ~  24^^  42  ~  22^^  44  ^  22^  42^^  23^^  34  “  44^^  32<^^  23  44^^  42^^  23^  34  24^^  42^^  22^^  44)  3  ■^'  ( 34^^  23<^^  44 

-  Cl  ^46/ 4 Vi  ^3  +  Cl  ’,4c/  42^i  :4<.i  33  ^  Cl  34a  22^^  33  (-i  34<^^  44  "  34<^^  22^^  32^'  23  “  C.I  x^il  22CI 42^^  24  “  24^^  32<^^  22^^  33  “  ^  24^^  32<^^  22^^  44 

-^i^4Ci3  3a  32^f44^^^  :4^i52^i23  +  ^/L^L32^02)^'4]^^^G2^Lu“^L3V^44)P22<  0  )  +  (  ^  22^/ 44  “  ^^2^^  24)  P 33^  ^  24^^32  “  22 ^ P44(  ^  ) . 

(A8) 

/^  =  .V  ~  ^  [ti  4^(  “  Cl  4Vi  ^  Vi  ^,4  “  a  -5  v’^  ^-.Ci  23  +  Cl  22^i  T,  "F  3^i  44 ^  3 V^  23^^ 44)  2  “  42^^  23^  22^^  33  “  32^^  23  22^^  44 ^  ^^3  44  42^^  24) <-  3  ^ 

^(  -aT/G2^^24^^^22^02^^23^Lu“^^^33^G2^^23^Lu”'^^24<^02‘^L'2^^23~^^22^^l3^  -<^^22^LV3^L32^i23“^G2t^i3)‘^‘4]”‘^O2^L33P:2(0) 

-f  Cl^2^l  2\P^^i  0  )  +  (  ii  22^^33”  ^L32^^  23)P44(  ^  )•  LA9) 

M  =  —  {cl22^  Cl  y.-T  ct^) .  (A  10) 

tV  =  ClizCl}}^  Cly2<^l2}'^  ClizCCu'^  Clx}Cl^4~  Cl2^Cl^2'  '  (All) 

L  =  a-,  <A12) 

where 

A  =  a  32i:i  23^^  44  ^  42^^  24^^  33  ”  42^^  23^^  34  ”  Cl  zz^i  \}Ci  44  • 
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